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ABSTRACT 


This  paper  presents  an  analysis  of  a non-linear,  single -axis  automatic 
flight  control  system  of  the  "flicker"  type.  The  investigation  was  initiated 
by  performing  theoretical  analysis  of  a simplified  system.  As  a result  of 
these  approximations,  it  was  possible  to  formulate  generalities  concerning 
steady  state  system  performance.  However,  this  method  is  limited  by  the 
difficulties  encountered  in  non-linear  analysis  and  therefore  electronic 
simulation  was  applied.  The  results  obtained  from  simulation  not  only  veri- 
fied the  findings  of  the  theoretical  aialysis,  but  also  extended  them  into 
more  difficult  areas. 

In  the  situation  considered  by  this  investigation,  the  controlled  body 
acts  in  a normal  linear  manner;  it  is  accelerated  proportionally  as  surfaces 
are  deflected,  and  it  may  be  represented  in  linear  transfer  function  terms. 
However,  the  "Flicker"  type  performance  of  the  automatic  control  system  is 
of  a non-linear  nature;  its  operation  is  preprogrammed  and  dependent  upon 
airframe  motion  only  for  selection  of  the  phase  of  the  program  to  be  followed. 
Thus  the  combination  of  control  system  and  controlled  object  is  non-linear 
in  its  overall  performance  even  though  most  of  its  components  are  linear  and 
may  be  represented  by  transfer  functions. 

A former  analysis,  AF  Memorandum  Report  MCREXG03,  dealt  with  the  steady 
state  frequency  response  of  a "flicker"  type  control  system  in  which  the 
actuator  traveled  at  constant  speed  while  changing  position.  The  inves- 
tigation presented  here  is  an  extension  of  that  study  into  a dual-mode  type 
surface  travel.  It  is  an  attempt  to  form  design  criteria  through  which 
system  simplification  and  improved  performance  may  be  obtained. 
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SECTION  I 


INTRODUCTION 


1.  The  purpose  of  the  following  investigation  was  to  determine: 

a.  the  roll  motion  frequency  response  of  a missile  equipped  with 
a Flicker-Type  Control  System  for  the  aileron  control  surfaces  hnd 

b.  the  significance  of  the  characteristics  and  the  mutual  relation- 
ships of  the  components  and  sub-comoonents  of  the  Flicker-Type  Control  System 
utilized  for  roll  stabilization  of  a missile. 

2.1  The  three  types  of  automatic  pilots,  commonly  used  for  aircraft 
control,  are  the  velocity,  displacement,  and  force  type  sjrstems  described  for 
example  in  pp  1S7-1&9  AF  Technical  Report  5677*  These  are  linear  control 
systems  in  their  fundamental  design,  i.e.,  the  corrective  control  exercised 
is  a linear  function  of  aircraft  motion.  These  systems  may  consist  of: 
sensing  elements  reacting  linearly  to  aircraft  motions  - for  example,  angular 
velocity  and  displacement  about  a particular  axis  - actuators  reacting  linearly 
to  sensing  element  outputs,  control  surfaces  reacting  linearly  to  actuators, 
and,  for  small  deviations,  the  airframe  reacting  linearly  to  control  surface 
deflections. 

2.2  In  contrast  to  the  commonly  known  linear  control  systems  as  men- 
tioned above,  a Flicker  Type  Control  System  is  one  in  which  the  servo  response 
is  not  dependent  upon  the  magnitude  of  the  command  signal  but  is  programmed 
according  to  the  sign  of  this  signal.  Figure  2.2.1shows  a generalized  block 
diagram  of  the  closed  loop  including  the  missile  and  a flicker  type  control 
system'  in  one  of  the  several  possible  arrangements.  The  plots  adjacent  to 

the  main  components  are  to  demonstrate  their  operational  performances  in  a 
qualitative  manner.  There  are  other  varieties  which  are  similar  to  this  type 
deviating  slightly  in  one  or  mere  sub-components,  for  exam-ole,  the  servo  res- 
ponse ray  depend  upon  the  amplitude  in  a stepwise  manner,  including  a zero 
bound  dead  zone. 

The  block  diagram  of  the  roll  control  system  considered  by  this 
report  is  shown  in  Figure  2.3.1.  The  roll  angular  velocity  of  the  motion  of 
the  missile  (1)  is  sensed  by  a gyroscope  (2)  and  its  rate  output  signal  is 
integrated  (3)  to  yield  the  missile  angular  displacement.  Both  values,  rate 
and  displacement  represented  as  electrical  signals,  are  added  in  a predeter- 
mined proportion  (4)  and  form,  the  "summing  command".  When  the  sign  of  the 
coi  bined  signals  (the  summing  command)  changes,  a discriminator  relay  (5) 
(sometimes  called  "phase  detector")  engages  the  actuator  (6)  to  drive  the  sur- 
fact  into  one  direction  with  a given  fixed  angular  velocity;  it  continues  to 
drive  the  surface  in  this  direction  until  either  a reversing  signal  is  received 
or  a physical  stop  is  encountered.  The  missile  reverses  its  direction  of  roll 
accordingly,  and  the  whole  process  repeats  itself  continuously.  For  this 
investigation  the  missile  is  considered  to  be  undamped  in  the  roll  axis. 
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In  a system  described  later,  the  surfaces  are  released  upon  change 
in  sign  of  the  mixer  signal,  driven  to  the  null  position  by  the  opposing  tor- 
que of  the  airstrear , and  then  motor  driven  to  the  opposite  side.  A system 
with  this  performance  is  called  the  "feathering  type"  flicker  control  system. 
Surface  deflections  introduce  proportional  torques  and,  since  the  missile  is 
assumed  to  be  undamped  about  the  longitudinal  axis,  they  can  be  considered  in 
terms  of  missile  angular  roll  acceleration.  Because  of  the  flip-flop  or  flicker 
characteristic  of  the  relay,  this  system  v;ill  have  a stead;/  state  of  oscillation. 

It  can  be  seen  from  Figure  2.3.1  that  there  are  two  major  sources  of 
non-linearities  in  this  type  of  control  system  which  do  not  generally  exist 
in  more  conventional  autopilots,  namely,  the  so  called  flip-flop  action  of  the 
relay  and,  because  of  certain  design  characteristics,  the  fact  that  the  control 
surfaces  ray  reach  the  physical  limits  in  their  travel.  This,  incidentally, 
is  a necessity  for  the  proper  performance  if  the  curve  of  the  surface  travel 
plotted  against  time  is  symmetrical  with  respect  to  a vertically  drawn  raddle 
line.  The  Flicker  Type  System  is,  therefore,  a basically  non-linear  type  of 
the  control  system. 

The  moment  at  which  the  mixer  signal  changes  sign  is  said  to  be 
the  "tripping  point"  or  "trigger  point"  i.e.,  the  point  at  which  the  control 
surfaces  should  begin  to  move  back  towards  the  neutral  position.  Assuming  a 
steady  state  of  the  missile  roll  oscillation,  the  surface  m.otion  is  a periodic 
(non-sinusoidal)  oscillation  whose  frequency  obviously  has  to  be  the  sam.e  as 
that  of  the  missile.  The  period  of  this  oscillation  is  determined  by  twice 
the  time  elapsed  between  two  consecutive  tripping  points.  The  period  of  os- 
cillation (for  a given  individual  system)  depends  on  the  phase  lag  of  the  rate 
gyro  output  signal,  relative  to  the  true  angular  roll  velocity  component,  and 
on  the  ratio  of  the  gain  factors  by  which  the  rate  gyro  signal  and  the  displace- 
rient  signal  (obtained  b-  integrating  the  rate  signal)  are  amplified  before 
being  fed  to  the  mixer.  Since  these  two  components  can  be  easily  varied  to 
gain  desired  reponses,  their  influence,  as  system  parameters,  is  Of  particular 
interest.  The  other  components  which  E:ay  be  adjusted  are  actuator  speed,  range 
of  surface  travel,  and  surface  effectiveness. 

It  will  be  seen  that  for  the  feathering  type  system  (when  the  curve 
representing  the  surface  travel  is  coimposed  of  straight  lines  and  when  physical 
stops  are  not  encountered)  the  amplitude  of  the  oscillation  is  affected  in  a 
linear  fashion  with  respect  to  these  latter  components  over  large  performance 
ranges,  but  that  these  have  no  effect  upon  the  frequency  of  missile  oscillation. 
However,  when  physical  stops  are  encountered,  both  the  range  of  travel  and  the 
actuator  speed  become  factors  in  determining  frequency.  The  absolute  magnitude 
of  the  summing  command  signal  (velocity  component  plus  position  component)  has 
no  significance  as  long  as  it  is  appreciably  greater  than  ths  threshold  value 
of  the  relay.  Hence  only  parameters  which  affect  the  change  of  sign  in  this 
mixer  signal  (summing  command),  in  respect  to  time,  are  to  be  considered. 
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SECTION  II 


REAC  SIMJLATION  DESCRIPTIONS 

3.1  It  is  well  known  that  an  attempt  to  perform  an  analysis  of  any 
system,  which  includes  pronouncedly  non-linear  functions,  is  necessarily 
tedious,  difficult,  and  time  consuming.  Certain  problems  involving  simple, 
periodic,  non-linear  functions  sometimes  lend  themselves  readily  to  a 
solution.  However,  the  non-linearities  involved  in  the  control  system 
under  consideration  do  not  fit  into  this  category  and,  hence,  mathematical 
analysis  of  this  problem  in  its  entirety  will  be  highly  complex;  the  final 
solution  would  very  likely  be  in  a form  which  would  render  applications 
difficult . 

However,  some  mathematic  solution  of  the  problem  is  desirable  as 
a guide  by  which  other  investigations  may  be  gauged.  Mathematical  analysis 
is  therefore  performed  for  a simplified  system,  and  yields  certain  valuable 
approximations  to  the  problem.  For  more  complete  and  quick  information  on 
missile  behavior,  an  electronic  analogue  computer  was  used  for  simulation. 
Most  of  the  complexities  of  the  system  were  included  in  this  solution.  The 
varying  of  parameters,  to  investigate  resulting  system  changes,  which  prove 
difficult  by  mathematical  means,  is  easily  accomplished  on  the  simulator. 
Hence  the  burden  of  this  investigation  was  borne  by  electronic  analogue 
simulation.  In  addition,  the  results  obtained  from  simulation  research 
were  frequently  a helpful  guide  in  the  attempt  to  find  generalized  mathe- 
matical formulations  or  graphs  for  approximate  solutions  of  problems  result- 
ing from  the  changing  of  parameters.  Such  graphs  can  be  used  advantageously 
for  design  work. 
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3.2  The  generalized  system,  as  it  is  considered  for  simulation  on 
a Reeves  Electronic  Analogue  computer,  has  been  represented  in  Figure 
2.3.1.  For  simulation,  the  missile  body  is  treated  in  the  manner  given  in 
Figure  3.2.1.  The  input  & representing  the  roll  acceleration  due  to  sur- 
face  deflection  is  integrated  twice  by  the  missile.  The  first  integration 
yields  the  missile  angular  roll  rate  component  to  which  the  rate  gyro 
(block  #2  in  Figure  2.3.1)  is  sensitive.  The  second  integration  by  the 
missile  yields  the  missile  angular  roll  displacement  which  is  not  directly 
utilized  for  control  purposes. 

The  characteristics  of  the  rate  gyro  (block  //2,  Figure  2.3.1)  are 
described  in  detail  by  Figure  3*2.2  (upper  part).  The  gain  of  the  gyro 
system,  i.e.,  volts  output  per  degree  sec*^-  input,  is  of  no  significance 
in  determining  the  zero  points  of  the  mixer  signal;  the  zeros  are  determined 
only  by  the  ratio  of  the  rate  signal  amplitude  to  the  displacement  signal 
amplitude,  as  will  be  shown.  However,  in  order  to  obtain  a comparison 
between  true  roll  angular  velocity  (input)  and  indicated  angular  roll 
velocity  (output)  for  evaluation  purposes,  this  gain  factor,  as  well  as 
certain  characteristics  of  the  simulator  setup  must  be  considered.  By 
assuming,  in  the  classical  representation  of  a gyro, 

O)  D + Bi)  * KD  - - g,  <f 

a steady  state  related  to  a step  input,  it  is  seen  that 

a)  2- -it? 

In  the  RHAC  setup  for  a rate  gyro,  seen  in  Figure  3.2.2,  (lower  part) 

c f*-abD  B=clc  K=  abd 

and  the  input  is  ; substituting  these  values  in 


VJADCTR  53-202 


2 


/*_  jsr  ' 

equation  (2)  yfelds,  (3)  f . Hence,  in  this  setup,  true  and 

indicated  angular  roll  velocities  are  related  by  the  factor  d/g,  ■ 
For  example,  consider  a gyro  which  has  a frequency  of  5.033  cps  and  a 
damping  ratio  ( ^ ) sr  0.31&J  then 

uh*  31,  62i  K*  **  * <*<>•  3*  Z1  “ 10 ■ 

It  can  be  seen  that  the  gains  of  the  setup,  which  determine  B and  K,  are 
the  individual  gain  factors  a and  c,  and  a,  b,  and  d respectively.  The 
operator  will  choose  a reasonable  combination  of  a,  b,  c,  d,  for  the  pur- 
pose of  simulation.  If,  for  instance, 

C = s a=  10  b=/o  c=l  d=<o ( thin  3=20  fa  1000 

< * 

and  the  factor  by  which  ^ is  to  be  multiplied  to  obtain  a steady  state 
output  signal  of  unity  magnitude  for  a step  input  signal  of  unity  magnitude 
is  d/G^. 


It  has  been  previously  stated  that  the  signals  utilized  for  control 

of  the  missile  are  "indicated  angular  roll  rate"  and  "indicated  roll  angle" 

obtained  from  the  farmer  by  an  integrating  motor  (block  #3,  Figure  2.3.1). 

The  sole  purpose  of  the  mixing  device  (block  #4,  Figure  2.3.1)  is  to  sum 

these  signals  continuously  in  a predetermined  proportion.  Schematic  details 

of  this  operation  are  shown  in  Figure  3.2.3  with  its  REAC  simulation.  Since 

only  the  relative  magnitude  of  these  signals  and  not  their  absolute  values 

affect  the  operation  of  the  control  system,  the  gain  of  the  rate  signal  is 

arbitrarily  considered  to  be  unity.  Thus  in  this  setup  the  ratio  of 

Gain  factor  of  the  rate  signal  is  simply  I/G20 

Gain  factor  of  the  displacement  signal 
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The  phase  detector  (block  #5,  Figure  2.3.1)  is  shown  with  two 
different  types  of  simulation  by  which  it  was  represented  in  Figure  3.2.4. 

In  spite  of  a comparatively  large  time  delay  inherent  in  the  computer , the 
REAC  simulation  for  the  discriminator  proved  satisfactory  for  the  non- 
feathering type  control  system;  it  was  not  satisfactory  for  the  faster 
operating  feathering  control  system  for  which  the  discriminator  was  simu- 
lated by  an  Armour  Special  Function  Generator  which  has  less  time  delay. 

The  output  of  the  discriminator  has  a very  good  rectangular  wave  form  with 

frequency  identical  to  that  of  the  mixer  signal  and  a constant  amplitude. 

• 

The  motor  actuator  (block  #6,  Figure  2.3.1)  drives  the  control 
surfaces  as  programmed  by  the  relay  output.  The  rotor  runs  at  constant 
speed,  and  hence  the  position  of  the  ailerons,  until  they  arrive  at  phy- 
sical stops,  is  the  integral  of  the  relay  si, [rial.  In  some  systems  (non- 
feathering type)  the  surfaces  will  contact  stops  in  normal,  steady  state 
oscillation;  in  others  (some  of  the  feathering  type)  they  may  contact 
stops  only  when  the  missile  is  subjected  to  disturbing  torques.  The  motor 
actuator  is  prevented  from,  operating  at  constant  speed  by  the  varying  tor- 
que of  the  air stream  and,  also  initially,  by  the  inertia  of  the  body. 

Thus,  while  the  deflection  of  the  aileron  is  considered  here  to  be  graphi- 
cally represented  by  a straight  line,  in  reality  it  will  be  slowed  slightly, 
as  deflection  increases,  by  the  increasing  opposing  torque  of  the  air- 
stream.  However,  because  the  actuator  output  is  large  with  respect  to  any 
load  placed  upon  it,  a general  picture  of  system:  operation  is  gotten  with- 
out considering  the  time  delays  in  the  actuator.  A later,  detailed  analysis 
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is  made  in  which  the  effect  of  including  the  time  delay  due  to  inertia  is 
investigated.  However,  wind  tunnel  tests  indicate  that  the  loading  effect 
caused  by  the  airstream  is  negligibly  small. 


The  RSAG  simulations  for  non-feathering  and  feathering  actuators 
appear  in  Figures  3*2.5  and  3*2.6  respectively.  The  following  paragraphs 
are  a detailed  description  of  their  operation: 

A.  Non-feathering  actuator:  The  REAC  setup  (Figure  3*2.5) 
simulating  a non-feathering  actuator,  with  no  time  lag,  is  an  integrating 
amplifier  connected  with  a limiting  device.  Their  combined  effect  is  as 
follows : 

a.  To  simulate  the  control  surfaces  being  driven  by  the  con- 
stant speed  motor;  the  relay  signal  is  integrated. 

b.  This  integration  continues  until  the  voltage  output  of  the 

f 

amplifier  is  equal  to  a preset  value  in  the  limiter  at  which  time  the  in- 
tegration stops  and  the  output  voltage  holds  constant,  simulating  the 
action  of  the  control  surface  as  it  contacts  a physical  stop. 

c.  IJhen  the  relay  signal  changes  sign,  the  integration 
process  resumes  and  the  output  voltage  of  the  amplifier,  which  has  been 
held  at  a preset  level,  changes  accordingly,  thus  simulating  the  surface 
being  driven  back  to  zero  and  deflected  to  the  opposite  side. 

B.  Feathering  System: 

A flicker-type  system  in  which  the  deflected  control  surface 
is  allowed  to  "free-wheel"  back  to  the  neutral  position  propelled  only  by 
the  opposing  torque  of  the  airstream’:  is  called  a Feathering  Flicker  Type 
System.  A simplified  analysis  of  the  feathering  system  is  made  by  consi- 
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dering  it  to  be  the  equivalent  of  a non-feathering  system  in  which  the 
outbound  actuator  speed  is  slower  than  the  inbound  speed. 

No  accurate  data  relating  to  the  aerodynamic  performance  of 
the  control  surfaces  when  free-wheeling  are  available.  Even  if  such  data 
were  available  the  coefficients  involved  would  assume  different  values  under 
various  conditions  such  as  different  airspeeds  and  densities;  also,  the 
clutch,  by  which  the  surface  is  engaged  to  aid  released  from  the  actuator, 
by  its  very  nature,  introduces  a variable  delay.  Hence,  because  of  the 
uncertainty  as  to  the  exact  dynamic  characteristics  of  these  surfaces 
during  the  feathering  period,  it  was  found  that  the  results  of  this  in- 
vestigation would  not  be  changed  significantly  if  this  second  order  sys- 
tem were  replaced  by  a first  order  system.  The  errors  introduced  by  such 
a simplification  are  not  greater  than  those  which  would  originate  in  se- 
lecting one  particular  fixed  quadratic  as  representative  of  the  entire 
range  of  surface  operation. 

A graphic  representation  of  the  "free-wheeling"  control  sur- 
face would  very  likely  show  it  to  be  a nearly  critically  damped  oscillation. 
The  spring  constant  supplied  by  the  airstream,  and  the  moment  of  inertia 
of  the  control  surfaces  are  such  that  the  surface  system  has  a high  natural 
frequency.  Thus,  the  time  taken  by  the  "free-wheeling"  surfaces  in  reaching 
the  neutral  position  when  returning  from  a deflection  is  comparatively  small 
and  independent  of  the  amplitude  of  deflection.  However,  this  time  deter- 
mines the  point  at  which  the  surfaces  are  again  engaged  to  the  actuator. 
Because  of  the  high  frequency  by  which  the  surface  system  feathers,  the 
integral  of  the  curve  representing  the  surface  travel  versus  time  during 
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the  feathering  time  will  be  relatively  small  compared  with  the  integral  of 
the  preceding  portion  covered  by  the  outbound  curvej  in  other  words,  there 
will  be  only  a small  change  of  the  roll  velocity  during  the  feathering 
operation.  Consequently  the  system  may  be  simulated,  without  sacrifice  to 
accuracy,  by  any  other  system  which  has  a relatively  small  integral  and 


would  reproduce  the  time  required  in  returning  to  the  null  position.  Con- 
sider, for  example,  an  exponential  of  the  form  ((*)=■  — 

where  0&  - surface  deflection  at  the  reversing  point,  a = constant,  and 

- / 

(L~  j y T = Time  constant.  Such  a function  may  be  made  to  be  zero  at  any 


reasonable  desired  time  for  a given  deflection.  The  amplitude  of  the  sur- 
face at  the  reversing  point  was  found  to  be  from  10°-  20°  for  most  con- 
ditions considered  by  this  investigation.  For  a fixed  value  of  the 

time  consumed  for  the  surface  travel  from  & & to  zero  will  differ  only  by 
a few  milliseconds  as  the  amount  of  deflection  under  different  flight  condi- 
tions varies.  This  variation,  however,  can  be  considered  negligible  for  the 
purposes  of  this  investigation. 

Until  the  relay  signal  changes  sign  the  simulation  of  the  feather- 
ing system  is  equal  or  similar  to  that  for  the  non-feathering  system.  Then, 
rather  than  returning  to  zero  along  a straight  line,  the  signal  returns 
along  an  exponential  path  and  after  that  follows  a straight  line  outbound . 
This  simulation  is  accomplished  in  the  manner  shown  in  Figure  3.2.6  which 
is  described  as  follows: 

The  output  of  the  Special  Function  Generator  "A"  representing  the 
discriminator  relay  (Block  #5  in  Figure  2.3.1)  has  been  previously  described 
to  be  dependent  for  sign  upon  the  summing  command  and  of  fixed  absolute 
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amplitude  C^j  . Assume  this  sign  to  be  positive;  it  is  fed  into  both 

the  limiting  integrating  amplifier  and  into  a summing  amplifier.  The 
integrating  amplifier  output  is  connected  to  a second  special  function 
generator  "B".  Its  output  has  the  same  amplitude  as  "A”  but 

its  sign  depends  on  the  integrating  amplifier  output.  The  output  from 
generator  "B"  is  also  fed  to  the  summing  amplifier.  Jf  , as 

assumed,  has  been  positive  for  some  time,  will  be  negative  and, 

as  their  amplitudes  are  equal  but  of  opposite  sign,  the  summing  ampli- 
fier outpi.it  is  zero  and  the  relay  controlling  the  integrating  amplifier 
feedback  remains  open.  Suppose  the  summing  command  changes  the  sign  and 
consequently  "A"  changes  its  sign.  Then  the  output  of  the  summing  ampli- 
fier becomes  suddenly  negative  and  closes  the  relay  in  the  integrating 
amplifier  feedback  thus,  according  to  the  above  equation,  causing  an  ex- 
ponential decay  of  the  integrating  amplifier  signal  which  passes  through 
the  null  position  to  a signal  of  opposite  polarity.  However,  the  very 
instant  when  the  integrating  amplifier  output  signal  changes  its  polarity 
the  output  of  generator  "B"  changes  its  sign,  too,  thus  producing  again  a 
zero  output  of  the  summing  amplifier  and  the  opening  of  the  feedback  cir- 
cuit allowing  the  integrating  amplifier  to  resume  its  integrating  function 
starting  from  zero.  Figure  3*2.7  is  a graphic  representation  of  the  signal 
sent  to  the  missile.  At  time  0 the  integrating  amplifier  begins  its  output 
at  Tq  generator  "A"  changes  sign  resulting  in  the  exponential  decay;  at  T2 
the  exponential  decay  has  reached  the  null  position  and  generator  "B" 
changes  sign  causing  the  relay  to  open;  at  T3  the  generator  "A"  changes 
sign  again,  etc. 
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Several  components  of  the  Flicker  Type  System  have  inherent 
time  delays.  Some  of  these,  even  though  they  are  relatively  small,  could 
be  significant.  Therefore,  the  RSAC  setup  shown  in  Figure  3*2.8  was  used 
to  investigate  this  possibility.  A desired  time  delay  could  be  introduced 
at  any  point  in  a setup  by  merely  inserting  this  arrangement  of  amplifiers 
into  the  circuitry.  The  setup  is  similar  to  that  normally  used  to  simulate 
a first  order  system.  Consider  the  time  delay  effect  for  example  on  an 
integrating  motor.  If  the  motor  is  assumed  to  have  negligible  moment  of 
inertia  so  that  it  may  be  represented  as  a pure  integrator  its  transfer 
function  may  be  expressed  in  standard  La  Place  terminology  as  1 , When  a 
step  voltage  is  applied  to  this  motor  the  result  is  such  that 
where  K is  the  magnitude  of  the  step;  in  the  time  domain  (B0  * K-t  . if 
the  moment  of  inertia  of  the  motor  is  considered  significant,  a time  delay 
results.  By  introducing  the  time  delay  the  transfer  function  of  the  motor 


becomes  * ^en  a step  is  applied  the  result  is  such  that 


or  in  the  ti 


time  domain:  ^0  ^ )j(  . 


Question  may  be  raised  concerning  the  factor  q in  the  setup  (i.e.,  the  q in 
the  numerator  of  the  transfer  function).  If  this  factor  was  omitted  the  re- 


would  be  @0  — jt  y,  (f  & ^ . It  can  be 


seen  that  such  an 


arrangement  would  introduce  a change  in  steady  state  value  of  the  motor  out- 
put speed.  Since  it  is  desired  only  to  intorduce  a time  delay,  the  "q  factor" 
is  accordingly  included  in  the  setup. 

The  Discriminator  Relay  has  been  mentioned  as  sensing  a change  in 
sign  with  ideal  accuracy.  In  reality,  a certain  magnitude  of  change  must 
occur  before  the  relay  becomes  cognizant  of  the  change.  This  magnitude. 
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necessary  to  obtain  a response  from  the  relay,  is  said  to  be  its  dead 
space.  A similar  dead  space  (due  to  play  in  the  actuator  linkage)  often 
occurs  in  mechanical  systems  possessing  a restoring  force.  This  dead  space 
was  simulated  by  the  REAC  setup  shown  in  Figure  3.2.9  with  Z0  being  the 
mixer  signal  as  it  is  sensed  by  the  relay.  The  change  of  sign  of  t1  is  sig- 
nal occurs  somewhat  later  than  the  true  change. 

SECTION  III 

THEORETICAL  CONSIDERATIONS 


4.1  A brief  comparison  between  the  performance  of  an  airframe  equipped 
with  a conventional  autopilot  and  a missile  controlled  by  a flicker  type 
control  system  will  help  to  clarify  the  problem  of  investigating  the  perfor- 
mance of  a flicker  type  system.  A conventional  autopilot  in  an  airframe 
reacts  to  a disturbance  in  the  manner  shown  in  Figure  4.1.1.  A disturbance 
occuring  at  t0  results  in  displacing  the  airframe  from  its  original  position. 
The  surface  is  deflected  instantaneously  by  the  autopilot  so  as  to  counteract 
this  deviation.  However,  before  the  airframe  reaches  its  zero  position,  t2> 
the  influence  of  the  velocity  signal  causes  the  control  surface  to  pass 
through  its  zero  position  thus  producing  a lead  of  control  surface  deflection 
relative  to  airframe  displacement.  The  airframe  overshoot  during  the  next 
half  cycle  is  then  greatly  reduced.  This  is  repeated  for  perhaps  one  more 
oscillation  (as  at  t^  and  t^)  until  both  vehicle  and  control  surface  are  again 
in  their  original  positions.  During  this  entire  process  the  action  of  the 
control  surface  was  continuously  dictated  by  outputs  of  sensing  elements  re- 
acting to  various  components  of  airframe  notion.  A vehicle  equipped  with  such 
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a control  system  is  said  to  have  a quiescent  zero  position. 

A system  consisting  of  a missile  with  a flicker  type  control  has 
no  such  quiescent  zero  position.  It  is  continuously  in  oscillation  even 
when  in  the  undisturbed  state.  The  magnitude  and  frequency  of  this  oscil- 
lation are  well  defined  and  maintained.  The  reason  for  this  behavior  is  that 
the  surface  motion  is  dependent  on  the  missile  motion  in  a discontinuous 
stepwise  manner.  It  responds  to  missile  travel  on  one  side  or  the  other  of 
a zero  position  by  surface  travel  to  an  extreme  position.  There  is  no  pro- 
vision for  the  surface  to  remain  in  a neutral  position  and,  hence,  even  if  the 
missile  were  to  assume  a zero  position  by  pure  accident,  it  would  be  forced 
from  this  state  by  the  surface  travel. 

4.2  From  what  has  been  said  previously  one  Eight  conclude  that  only  a 
signal  which  is  a linear  function  of  missile  position  is  necessary  for  a 
flicker  type  control.  However,  a system  using  only  an  angular  position  sig- 
nal is  shown  by  the  following  consideration  to  be  unworkable.  Consider  a 
system  in  which  a sensing  device  detects  change  of  missile  position  to  either 
side  of  a zero  position.  This  activates  an  actuator  to  drive  the  control  sur- 
face at  a constant  speed  to  the  extreme  position  which  tends  to  return  the 
missile  to  the  opposite  side  of  zero.  Assume  that  such  a system,  be  workable. 
It  would  result  in  the  periodic  acceleration  curve  shown  in  Figure  4*2.1. 

Such  an  acceleration  curve  would  yield  the  position  curve  also  shown  in 
Figure  4*2.1.  The  operating  condition,  i.e.,  the  position  = zero  at  tripping, 
is  shown  to  be  feasible  only  if  the  surface  travels  from  one  extreme  to  the 
other  in  zero  time  as  shown  in  Figure  4.2.2.  This  is  practically  unworkable 
and  is  rendered  more  so  by  the  fact  that  any  detector  used  would  contribute 
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a phase  lag  to  the  position  signal. 

Assume  next  a system  which  triggers  on  zeros  of  velocity  only. 

Such  a system  if  workable  rust  present  a periodic  acceleration  which  would 
yield  the  velocity  curve  shown  in  Figure  4«2.3»  It  can  be  seen  that  the 
velocity  zero  must  occur  in  the  middle  of  the  surface  dwell  period.  Hence, 
such  a system  ca-’  work  only  if  the  trigger  line,  i.e.,  zero  point  of  velocity, 
coincides  with  the  midpoint  of  the  surface  dwell;  this  is  possible  only  if 
surface  dwell  time  is  zero,  as  shown  in  Figure  4.2.4.  However,  such  a system 
does  not  sense  minor  deviations  which  occur  at  a very  low  rate,  and  the 
system  will  drift  from  the  zero  line. 

Since  angular  velocity  and  position  are  the  most  reliably  obtained 
components  of  missile  motion,  one  next  considers  a possibility  of  using  a 
combination  of  roll  angular  velocity  and  angular  position  to  trigger  the 
system..  Consider  a system  which  triggers  whenever: 

Gain  factor  (g-^)  x Velocity  Signal  + Gain 
Factor  (g2)  x Position  Signal  = 0. 

If  such  a system  were  to  operate  ideally,  the  components  of  missile  motion 
would  have  the  relationship  shown  in  Figure  4*2.5.  The  operation  of  this 
system  is  not  as  immediately  obvious  as  were  the  two  preceding  systems  but 
can  be.  ascertained  only  by  analysis.  Hence,  a theorectical  analysis  of  this 
system,  is  presented  in  paragraph  4.3. 

4.3  The  following  is  a theoretical  analysis  of  a somewhat  generalized 
flicker  type  control  system.  The  system,  considered  operates  as  follows:  The 
control  surface  is  driven  out  from:  the  neutral  position  at  K-|_  degrees  per 
second;  it  continues  outbound  until  (l)  it  reaches  physical  limits  or  (2) 
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it  reaches  a position  at  which  the  velocity  component  is  some  predetermined 
fraction  of  and  opposite  in  sign  to  the  position  component.  In  either  case 
when  gain  terms,  in  a ratio  inverse  to  this  fraction,  are  applied  to  the 
angular  velocity  and  position  signals,  and  when  the  sum  of  the  resulting 
signals  equals  zero,  the  system  triggers.  At  this  moment  the  control  surface 
returns  to  neutral  at  K ^ degrees  per  second,  is  driven  to  the  opposite  side 
at  K-j_  degrees  per  second,  and  the  cycle  repeats.  The  following  nomenclature 
is  used  for  this  analysis: 

- Outbound  surface  velocity.  (Deg/Sec) 

K2  r Inbound  surface  velocity.  (Deg/Sec) 
gq  r Gain  term  on  the  velocity  signal.  (1) 
g2  = Gain  term  on  the  position  signal.  (l/Sec) 

R * £2/ sl*  [*?<.] 

z - Time  to  drive  surfaces  from  neutral  to  stops.  (Sec) 
tit  r Time  surface  spends  at  stops  in  each  half  cycle.  (Sec) 

P aKj/fcj.  (1) 

V = Kissile  Angular  Velocity.  (Deg/Sec) 

V*  - Indicated  Kissile  Angular  Velocity  (Deg/Sec) 

D = Kissile  Angular  Position  (Deg) 

D*  = Indicated  Kissile  Angular  Position  (Deg) 

Vmax  = Maximum  velocity.  (Deg/5ec) 

Djaax  m Maximum  displacement  (Deg) 

V rp  - Velocity  when  system  triggers  (Deg/Sec) 

D/y*  s Position  when  system  triggers  (Deg) 
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t = Time  at  wh:  ch  velocity  is  zero.  (Sec) 

T s Time  at  which  tripping  occurs  (Sec) 

This  analysis  is  based  upon  the  following  assumptions: 

1.  The  system  is  in  steady  state  oscillation. 

2.  The  missile  is  assumed  to  have  no  damping  about  the  roll  axis. 

3.  At  tripping:  V,j  = RD-j.  or  g^V<p  ■ S^T  • 

4.  ^ K,  i.e.  PSj 

5.  Surface  effectiveness  is  assumed  to  be  l°/seca’per  degree  surface  deflection. 
Effectiveness  may  be  introduced  into  these  considerations  by  multiplying  the 
surface  velocities  Kq_  and  K2  by  the  effectiveness  factor.  Missile  response  will 
be  investigated  with  respect  to  missile  frequency  and  amplitude  of  oscillation 
as  they  are  related  to  P,  R,  V,  and  'C  • 

Surface  deflections,  in  the  system  described  above,  will  impart  a 
torque  to  the  missile  which  is  proportional  to  the  amount  of  deflection;  the 
missile  roll  acceleration  is  proportional  to  this  torque.  Disregarding  time 
delays,  a single  cycle  of  roll  acceleration  may  be  represented  as  shown  in 
Figure  4.3.1.  Roll  acceleration  is  obtained  by  multiplying  surface  deflection 
by  effectiveness.  Symbolic  representation  of  a cycle  of  surface  deflection  is 
given  in  the  following  six  equations: 

A = if  ost^'c 

k = If  r *t*0+*)x 

3.  A-K  r-tt  ;f  (•*«**)* 

4.  if  (i+H+P)T  * t (i+ M+-PH 
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5.  k-  - K,  v if  U+M+?jrsts  (z<-zm-?)v 

6.  if  (Z  + lMi-V)T  S is  l(l+M+p)T 

The  time  intervals  listed  above  are  hereafter  designated  interval  1,  2,  3* 
etc.  in  the  order  listed:  t]_  = 0 at  the  beginning  of  interval  1,  t2  **  0 

at  the  beginning  of  interval  2,  etc.  "t"  represents  continuous  time  from 
the  beginning  of  the  cycle. 

The  roll  velocity  is  the  integral  of  roll  acceleration  which  is 
symmetrical  about  the  time  axis  (A  e 0).  Hence,  the  roll  velocity  will  be 
symmetrical  about  an  axis  parallel  to  the  time  axis  and  since  the  missile 
is  assumed  to  be  controlled  about  a given  point,  this  axis  must  be  the  line 
V = 0.  Roll  velocity  is  described  by  the  following  equations: 

7*  V jKt,dt,  + C - - Ktar  >f  osi*  v 

(Where  it  is  assumed  that  A(o)  = 0 and  hence  V(o)  - VDlax  or  Vniin  and  is  ar- 
bitrarily taken  as  Vm^n  which  by  symmetry  is  - Vffiax) 

a.  V*  K v tx  + _ VthM  r<^  ct+N)t 

9.  (/=  t3  - Qfl  + + K,  MtJ-  I/**,  if  (Mr  *t* 

(t+  rn-?)r 

10-  V=  k*.«,  - if  (J+M+v) r < ts  (2+M+  pjt 

l/=  Kwcv  - ~fr  — Kitts  if  (2+M+rjt ■s.t s Cx+2M+p)T 
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la.  (/=  ^ 

if  il-tXb \+?)r*t'S  (Z+ZM*-ZP)t 


When  t = ((+ti+v)zf  l £«*  (See  Figure  4.3.1)  which  when  substituted  in 

(9)  yields: 


14 4,  = i K 1 V KM  rV  { K,  Pt 1 - V*,*v 


from  which:  (13)  . / 

*th  ay 


i K,tl(i+iM+p) 


Also  from  Figure  4.3«1>  it  can  be  seen  that  tripping  occurs  at  t = (t+M)Z 
which  when  substituted  in  (8)  yields: 

Mp  ~ ^ ft  V* «.x 

from  which:  (14) 

Vr  ■=  i(Ktl)U+in~?) 

tQ  is  defined  as  the  point  at  which  V = 0.  This  is  the  point  at 
which  the  area  under  the  acceleration  curve,  from  t = 0 to  t » (l+M+V) Z 
is  divided  in  half.  Examination  of  the  curve  and  considering  a previous 
assumption  that  K2  i Kp  reveals  that  t©  could  lie  in  either  the  first  or 
second  time  interval. 

Case  I If  t0  is  in  the  first  time  interval,  it  may  be  substituted  for 
tx  in  (7).  This  yields:  q ^ £ K,t,  - OK 

(IS)  t0  = rf[(j7zM+p) 

Case  II  If  tQ  is  in  the  second  tire  interval  it  can  be  substituted  in 
(8)  for  tg  • This  yields: 

K^tc-  Kay 
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or:  (16)  where  "t0"  is  calculated  from  the  beginning 

of  the  second  time  interval  so  that  the  actual  time  is^a  £ o + • 

The  Displacement  is  the  Intergral  of  Velocity  with  its  minimum  and 
maximum  at  zero  velocity.  Therefore  the  greatest  change  in  displacement 
Dmax  “ ^min  represented  by:  C 

= J 

bo 

Velocity,  a function  symmetrical  about  V * 0 and  since  the  missile  is  assumed 


to  be  in  a steady  state  of  undisturbed  oscillation  it  may  be  assumed  that  the 
point  about  which  it  oscillates  is  D = 0 and  that  displacement  will  be  sym- 
metrical about  this  point.  Therefore:  . 

0+Mt?)Z  +te 

1?-  JW  = - I J VCiJdt 

Displacement  at  the  tripping  point  will  therefore  be: 


18.  For  Case  I:  (/+M+ PJZ  +ta  t f / U 

2>t  =-jr  ■ ( Vdj&i-J  Jl/(tjdt 4 

19.  For  Case  lit  *'  , ,, 

j>T  = -%  J vctjdt  + J V(ydtx 

to  to 

These  are  shown  to  be  identical  by  the  following  considerations: 

t -HZ 

Expanding  equation  (18):  ~T)~«  » •—  2 f \Z(tt)dtf  Jtf  V 

7>*  i **  * ° pz 

- i JvNA  - if  WA  + / 


19.  For  Case  II: 


Gathering  terns  and  using  the  fact  that  VCtJ  - - V(tj 

r **  fT 

20-  j)  = f JWJA  +ij MJA  -i J vpMs 
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Expanding  (19)  in  * similar  manner  — ~~2  j VCtj)dt^ 

. ' ■fl'(t4)Jt4  - if vcts)dtc  + J hw&i 

*'Ja  o to 

and  using  the  fact  that  Vdj)  — "~VCty  and  V(t:)  = - veto 

x 

(21)  dr-ijvM 

o o o 

(20)  * (21)  i.  J)ft  * J)Ti 

Evaluating  for  J)t  by  substituting  values  for  mj  etc.  in  equation  (20) 
XT  Afr  2 

j)r  - QA+tfvO1?  - H.W 

-i j(«  **,  KHz2  - U4 

-J\  _ 0 V»ng  Z(hti-T)  . K MxJ  , K,Pzl  _ K.  MPz1 

JJf » ~ - j ^ 4 4 l 

J)r  • -n-p-2nr-{-n‘r%] 


(M)  J)a 


K,r> 


a+s+yp+T+V') 


The  introductory  description  of  the  system  contains  the  condition 
for  tripping  = -aa.  . Substituting  from  equation  (14)  and  (22): 

tei(t+ZM-p)~  ’Zfj£-(l+2t1+9MP+(,r+T>1) 

(23)  From  wM.k  « - &+**<■*>*  I'f-t+iP  .„a  „ 6(I*IM~P) 

a-  9fi*P-3 Rr  'RO+3SH  9MP+6P+  ?*) 
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It  is  of  considerable  interest  to  note  here  that,  idiile  surface 
speeds  were  significant  in  calculating  actual  values  of  velocity  and  dis- 
placement, they  enter  into  frequency  considerations  only  insofar  as  they  deter- 
mine the  time  taken  in  reaching  physical  stops  (t).  Two  cases  ef  particular 
interest  are  the  case  involving  no  physical  stops,  i.  e.,  M s d , and  the  case 
where  outbound  surface  velocity  equals  inbound  surface  velocity  i.e., 

K,-Ki  ~ 7W 

Case  where  H:0  (no  physical  stops) 

^ 6 0 - SL 

(24)  Z **  IKl+bp+V*) 

Here  it  .can  be  seen  that  if  H « 0 and  P « 1 the  system  cannot  operate.  Case 
where  P ■ 1 (one  inbound  and  outbound  surface  velocity) 

(25)  3(7-  fij  0r  ^ 

The  period  of  oscillation  in  the  general  case  is,  from  inspection  of  Figure 

4.3.1,  20+M  + P)T  • For  the  special  cases:  if  P = 1 the  period  is 

Ux+m)t  and  if  M - 0 the  period  - Z0+P)T.  The  maximum  oscillation 

of  the  missile  is  shown  in  (17)  to  be:  _ . , 

(l-Hi+P)Z  + ta 

Dtoctx  ® ~ T J V(t)  dt 

tc 

where  t0  is  the  time  at  which  the  angular  roll  velocity  is  zero.  As  has  been 
discussed  earlier  in  this  analysis  t0  = ti  or  t2»  tQ  has  been  calculated  in 


equation  (15)  and  (16).  In  Case  I,  if  t0  r % , t0  ■ Tr  \/ 

But  since  t0  is  not  easily  determined  it  is  desirable  to  state  these  conditions 
without  reference  to  t0  • Using  the  condition  for  Case  I that  tQ  - £T  we  have 


26.  £ > 


1 + X M + T> 

~ 


*•  / 2 2M+P 


Similarly  for  Case  II,  tc  * t^  , and  from  equation  (16)  tQ  « (- U2.tU-V)lf 
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To  have  any  meaning  £ 0 and  hence: 

27.  / Ss  2M+P 

Calculating  the  maximum  amplitude  for  Case  I 

n .-/>  . ,jl. 

J-'hULlt  J'tH" 1 ^ / 


Jmdt 


% Mv  rc  ^ 

J>h,ay  =i  J to  +ijmdts  tijmjJt* 

{.  0 0 0 

_ K^riiPni+xrn-pf-i-iHi-inP+P] 


28  • Aa*  = -^T 


Likewise  determining  the  maximum  amplitude  for  Case  II  in  which  , 

!<2M  + P andt 0.  (-l  + M + PJf  . n- 

1 )«av  * X ^ J'vfijdt  = if V(tt)dti  + 

^ to  t0  to  0 

+ {(  v(t+)dt+  + {JV(ts)dts 

0 o 

il.  A„x  - -frY 7+24M+  HP*  M MPi-  7 P*+ 12  Ml) 

Now  consider  the  system  to  have  a steady  state  in  which  the  sur- 
face dwells  in  its  extreme  position  for  part  of  the  cycle,  i.e.,  M 3t  0 # 
If  we  consider  the  two  extreme  values  for  P,  P B 1 or  P B 0,  we  get  the 


following  equalities: 


30.  From  equation  (25)  if  P = 1 M — 

/Pc  ""  () 

31.  And  substituting  P = 0 in  23.  » 2(4~  fa) 

in  30.  M y 0 only  if  0 < Rz  ^ / 
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Consider  £ to  be  some  fixed  time  as  it  will  be  in  the  practical  case; 

then,  as  R varies  from  0 to  p , M increases  in  value  and  the  frequency 

decreases.  In  31*  H>0only  if  4,<T$Z  —6  or  £ < 7?  ~ — • 

V 

Again  considering  V fixed  it  is  seen  that  varying  R from  ^ to  ^ results 
in  decreasing  M and  increasing  the  frequency.  If  one  considers  that  the 
effect  of  varying  R is  exactly  reversed  for  extreme  values  of  P,  it  occurs 
that  there  must  be  some  sort  of  a crossover  value  of  P.  Considering 
T — the  values  of  R which  yield  /V 13  oo  or"  0 c.fis.  , and  M = 0 

or  no  dwell  time,  are  plotted  against  varying  values  of  P in  Figure  4*3»2. 

It  is  seen  that  the  "infinite  dwell"  line  and  "no  dwell"  line  intersect  at 
about  P = 0.13.  This  is  superficially  interpreted  to  indicate  a condition  in 
which , if  the  surface  does  contact  its  stop  prior  to  reaching  a triggering 
condition,  the  system  will  not  operate.  Granting  that  this  value  of  P may  be 
considerably  altered  for  practical  cases,  it  can  be  assumed  that  the  point 
still  exists  and  will  be  a critical  point  in  system  operation.  It  is  noted 
that  for  the  value  'V  ~ this  point  occurs  very  near  to  the  value  of  P 

which  could  be  yielded  by  a feathering  system j the  symmetrical  case,  i.e., 

P as  1 is  far  away  from  this  point.  Furthermore,  by  differentiating  in 
(23) 

VT?  ~ VX-W'cP-l'R'c)1 

from  which  it  can  be  seen  that,  while  varying  X does  vary  the  magnitude  of 

^3  AT 

/j>  for  any  given  value  of  R,  varying  X does  not  change  the  value  of  P 
for  which  ^ 0 i.e.,  the  changeover  point.  Since  the  investigation 

of  the  subject  under  discussion  has  not  been  carried  any  farther  in  an  exten- 
sive manner,  it  should  be  kept  in  mind  that  the  picture  may  change  signifi- 
cantly when  practical  limitations  are  considered.  However,  there  is  a certain 
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indication  that  even  in  the  practical  case  such  a point  might  exist  and  should 
be  approached  with  caution. 


Reference  to  equation  (24)  shows  that,  for  the  "no  dwell"  condition, 
decreases  as  R increases.  That  is:  increasing  R increases  frequency. 

It  is  noted  that  the  condition  placed  upon  the  overall  system,  Vr  - - 1?J)T 
eliminates  the  sjnrmetrical  case  with  no  dwell. 

The  following  are  the  more  important  equations  developed  by  these 
considerations  where  the  surface  effectiveness  is  l°/sec2  per  degree  surface 
deflection.  Only  equations  involving  Kq  will  be  affected  by  changing  the 
effectiveness.  This  would  be  accomplished  by  substituting  C for  where 
C s EKj  j E being  the  effectivity  factor. 


V**,  - i K,  t*  ( ^ 2n  + p) 

Dtoix  = -§ 'f*  3MP+  P*J  if  tZZH+r 

D*,ay  = (7+  Z4M  + UP  + MMP+7PirUMt)  ift*2n+P 

I ff  *=  j 1C,  x 3 (!+  2M-P) 

Dt  = -*&(t+3PH-9MP+6P+Pl) 


V = 

M = 
M = 

r = 


o 


vh<«  t- 

UI,IH  t * (3+2M+P)  f- 
Pr  ± 

a-  <)VtP  -3Pv 


if  t>ZM+P 
if  t - ZM+V 

y titer  ally. 


Vr  -6 

m-M 

da-?) 

1?C<i-6P+?Z) 


if  p=0  M 
if  M~0 


l£t 

3(/-P t) 


Nett:  P-0  — ideal  feathering.  P«/  — Symmetrical  case./ 
M-0  S ho  physical  stofs. 
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By  the  method  of  approach  used  here,  it  appears  that,  for  the 
feathering  case,  the  frequency  increases  as  the  gain  ratio  is  increased. 

This  is  shown  by  the  formula  (24)  (when  M = 0)  but  it  is  more  clearly  seen 
by  reference  to  Figures  4.3.3,  4.3.4>  and  4.3.5  which  show  an  approximation 
of  the  Feathering  Case  at  0.5  cps,  1 cps  and  2 cps.  These  calculations 
assumed  ideal  sensing.  For  the  symmetrical  case  this  assumption  ignored  only 
a phase  lag  introduced  by  the  sensing  device;  this  will  be  readily  remedied 
in  paragraph  4.4.  However,  in  the  feathering  case  practical  sensing  devices 
will  not  merely  contribute  a phase  lag  to  the  velocity  signal  but  will  also 
modify  the  wave  form  significantly.  Hence,  one  must  be  careful  in  drawing 
conclusions  from  this  analysis  concerning  the  feathering  case.  It  is  believed 
that  good  results  could  be  obtained  by  analyzing  this  case  by  Fourier  methods. 
Once  an  approximation  has  been  made  of  the  acceleration  curve  by  a Fourier 
series,  the  velocity  and  displacement  series  are  gotten  through  termwise  inte- 
gration of  the  original  series.  The  fundamental  term  and  all  significant 
harmonic  terms  of  the  angular  velocity  are  subjected  to  the  dynamics  of  a rate 
gyro,  and  the  rate  gyro  output  is  obtained  as  a series.  This,  when  integrated, 
yields  the  indicated  position  signal.  These  ixiicated  components  of  missile 
motion  are  the  controlling  factors  in  missile  oscillation;  just  as  a group  of 
curves  are  assembled  later  in  the  report  for  the  VPR  chart,  in  Figure  4.4.4, 
curves  representing  several  frequencies  for  the  feathering  qystem  could  be 
assembled.  The  focal  point  of  such  a drawing  would  again  be  the  triggering 
line,  which  could  be  moved  as  described  in  4.4  to  represent  other  phase  lags  in 
the  system.  In  this  case  however,  the  gyro  used  would  be  an  integral  part  of 
the  chart  which  for  the  feathering  system  would  differ  in  uniformity  of  appear- 
ance from  the  VPR  chart,  but  should  be  at  least  as  useful  for  design  purposes. 
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4.4  The  analysis  in  4.3  shows  the  final  system  proposed  in  paragraph 
4.2  to  be  workable  within  the  stated  operating  limitation.  This  analysis 
assumed  subcomponents  without  time  lags.  Practically,  however,  such  a system 
is  not  feasible.  If  one  is  to  establish  design  criteria  on  this  system  a 
Bieans  must  be  devised  for  including  these  phase  lags. 

In  order  to  obtain  a clearer  picture  of  the  system  analyzed  in 
paragraph  4.3  by  equating  P g 1.,  graphic  representation  is  made  of  the  ideal 
system  action  at  three  operating  frequencies.  In  Figures  4.4.1,  4.4.2,  and 
4.4.3  the  angular  acceleration,  angular  velocity,  and  angular  displacement 
imparted  by  a so-called  symmetrical  flicker  type  system,  are  shown  in  their 
relationship  to  each  other.  The  figures  represent  frequencies  of  1.25  c.p.s., 
1 c.p.s.,  and  0.833  c.p.s.  respectively. 

These  curves  were  drawn  from  theoretical  considerations  presented 
in  paragraph  4.3.  Surface  stops  were  arbitrarily  set  at  t 15°  with  surface 
velocity  taken  as  120°/sec.  Evaluation  of  the  angular  velocity  and  displace- 
ment at  the  time  the  system  is  triggered  reveals  ratios  of  displacement  to 
velocity  to  be  1:5,  1:4.32,  and  1:333  for  frequencies  0.833#  1,  and  1.25 
respectively.  Hence  the  ratios  of  the  gain  applied  to  these  signals  as  they 
enter  the  mixer  must  be  the  inverse  of  this  to  obtain  triggering  at  this 
point.  The  broken  curve  represents  a comparison  of  the  rate  and  displacement 
signals  entering  the  mixer  if  these  gain  ratios  are  applied.  It  can  be  seen 
that  as  the  "gain  in  displacement"  is  increased  for  a given  gain  in  velocity, 
the  frequency  decreases. 

Any  component  which  tends  to  delay  these  signals  simultaneously, 
without  distorting  them  significantly  (for  example  the  rate  gyro),  can  be 
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represented  by  shifting  the  "trigger  line".,  on  these  graphs  a comparable 
distance  to  the  left.  A component  which  delays  one  signal  can  be  shown  by 
shifting  the  trigger  line  only  for  that  signal;  if  this  is  accompanied  by 
an  attenuation  the  graph  signal  must  also  be  attenuated  accordingly.  These 
statements  are  meant  to  be  indicative  of  approximations  which  can  be  made; 
they  ignore  filtering  actions  of  the  rate  gyro  which  change  the  shape  of  the 
velocity  curve  and,  hence,  of  the  position  curve.  A time  delay  in  the  inte- 
grating motor  component  makes  the  position  curve  even  more  sinusoidal.  For 
the  symmetrical  ease,  the  curves  are  already  very  close  to  a sinusoid,  thus, 
(apart  from  the  shift  due  to  phase  lag)  the  error  made  by  utilizing  the 
original  velocity  and  position  curves  instead  of  corrected  ones  is,  for  the 
purposes  under  consideration,  negligibly  small.  These  methods  should  prove 
very  advantageous  in  system  design  problems. 

The  only  truly  significant  parts  of  Figures  4*4.1,  4.4.2,  and  4.4.3 
are  those  quarter  cycles  immediately  preceeding  the  triggering  line*  A graph 
comparing  angular  velocity  and  positions  in  this  region  for  the  several 
listed  frequencies  is  shown  in  Figure  4.4.4.  For  these  curves,  aileron  stops 
were  placed  at  ±20°.  The  acceleration  curve  is  omitted  because  it  is  not 
directly  significant  to  these  considerations.  The  scales  for  velocity  and 
displacement  are  different  to  facilitate  closer  approximations.  Since  the 
acceleration  curve  preceeding  the  trigger  line  is  of  constant  value,  i.e.; 

20°/ sec^  (for  the  surface  is  dwelling  at  its  stop),  the  velocity  curve  from 
its  zero  point  to  the  trigger  line  is  a straight  line  with  a slope  20°/sec/sec. 
From  the  trigger  line  on  the  acceleration  curve  is  a straight  line  with  a slope 
of  -120°.  sec  ^/sec.  Therefore  from  this  point  to  its  maximum  the  velocity 
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curve  is  parabolic,  increasing  over  its  trigger  position  by  1.67  degrees. 

The  zero  point  of  the  displacement  occurs  at  the  point  of  maximum  velocity 
and  the  position  values  at  maximum  and  triggering  points  may  be  obtained 
by  reference  to  formulas  No.  29  and  22  in  paragraph  4.3. 

The  major  sources  of  phase  lag  and  attenuation  are  the  following: 

1.  Rate  gyro  which  results  in  delaying  and  amplifying  both  signals. 

2.  Integrating  motor  which  delays  and  attenuates  and  position 
signal  only. 

3.  Time  lag  in  the  relay  motor  - time  taken  by  the  relay  signal 
to  change  from  one  extreme  to  that  of  opposite  polarity.  This 
results  in  delaying  actuator  action  and  thus  causes  triggering 
to  occur  after  the  time  at  which  triggering  conditions  exist 
i.e.,  after  Vrp 

4.  Dead  space  in  relay  sensing  - results  in  a delay  in  beginning 
the  reverse  of  relay  signal  and  hence  also  delays  the  actual 
triggering. 

5.  Time  lag  in  the  actuator  - results  in  a delay  of  the  trigper 
line,  in  the  same  manner  as  3 and  4 above. 

The  only  one  of  these  delays  which  deals  with  one  signal  by  itself  is  that  in 
the  integrating  motor.  The  amount  of  phase  lag  and  attenuation  of  the  position 
signal  caused  by  the  integrating  motor  is  frequency  dependent  as  indicated 
by  Figure  4. 4.4.1. 

The  actual  sensing  device  producing  the  control  signal  for  the 
missile  ailerons  is  a rate  gyro  yielding  the  indicated  angular  roll  velocity 
signal.  This  signal  serves  a dual  purpose;  it  is  first  fed  into  the  summing 
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or  mixing  device  and  secondly  it  controls  an  integrating  motor  whose  output 
is  also  fed  into  the  summing  device.  Computer  study  reveals  that  the  intro- 
duction of  relatively  large  time  de]ys  into  the  integrating  motor  component 
has  a comparatively  small  effect  on  the  frequency  of  missile  oscillation.  In 
certain  instances  improvements  in  system  performance  result  from  the  intro- 
duction of  a time  delay.  These  results  are  rather  surprising  and  demand  an 
explanation. 

For  the  following  considerations  the  indicated  angular  roll  velocity 
signal  is  the  basic  signal.  Furthermore,  the  velocity  is  assumed  to  have  a 
true  sinusoidal  nature.  The  Figure  4.4.5  is  drawn  in  such  a way  that  the 
velocity  component  is  represented  by  a vertical  vector  V lying  in  the  positive 
imaginary  axis.  Then,  an  ideal  integrating  motor  component  with  no  time  delay 
would  yield  a signal  represented  by  a horizontal  vector  P lying  in  the  positive 
real  axis.  To  simplify  the  representation  of  the  components  relationship,  this 
vector  shall  be  of  unity  length.  For  a given  frequency  ( tt  in  radians/sec)  by 
which  the  missile  oscillates,  the  vector  V then  has  the  length  of  units  if 
the  gain  ratio  of  velocity  component  over  position  component  is  l:k.  (The 
diagram  in  Figure  4.4.5  is  drawn  for  a ratio  1:2  andfc*  = 8.)  The  resultant 
vector  R represents  a summing  command,  whose  zeros  determine  the  trigger  points 
of  the  control  system.  The  "time  line"  t-t  rotates  in  a clockwise  direction 
about  the  origin  of  the  complex  coordinated  0.  If  the  position  of  this  time 
line  is  chosen  such  that  it  is  perpendicular  to  the  resultant  vector  R,  the 
projections  of  the  vectors  V and  P on  this  time  line  are  of  equal  length  and 
opposite  in  sign.  This  is  the  condition  necessary  for  triggering  the  system. 

The  angle  between  t-t  and  OP  which  is  the  same  as  between  OR  and  0V  is  the 
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phase  angle  between  the  zero  of  the  velocity  component  and  the  instant  when 
triggering  occurs.  Consider  now  the  case  in  which  the  iriegrating  motor  has 
a noticeable  moment  of  inertia  represented  by  a certain  time  constant  in  its 
response  transfer  function.  The  integrating  motor  output  signal  then  lags 
behind  the  true  integral  of  the  velocity  component.  If  the  frequency  of 
missile  oscillation  is  assumed  not  to  change  appreciably,  the  locus  of  the 
position  vector,  as  the  time  constant  varies  from  zero  to  infinite,  is  very 

nearly  the  semicircle  PP-jPgO.  The  phase  lag  angle  is  given  as  tan  -jrf  • 

* 

where  T is  the  time  constant  of  the  motor  integrator.  As  can  be  seen 
by  graphical  representation,  for  a relatively  large  change  in  q there  is 
only  a very  slight  change  in  the  angle  VOR.  The  locus  of  the  vector  R,  as 
q varies  from  00  to  0,  is  the  semicircle  RR]_R2V.  Expressed  in  a qualitative 
manner:  as  to  the  response  of  the  system,  for  an  increasing  Time  Constant  the 
increase  of  phase  lag  (between  the  time  integral  of  the  indicated  angular 
velocity  component  and  the  actual  delayed  position  sigial)  will  be  cancelled 
by  the  attenuation  of  the  position  vector  thus  resulting  in  only  a small 
change  of  the  overall  missile  performance.  If,  for  example,  the  time  delay 
in  the  rotor  integrator  is  such  that  the  position  vector  0P]_  lies  in  the 
direction  of  the  time  line  t-t,  as  drawn  in  the  figure,  it  can  be  seen  that 
the  resultant  vector  has  exactly  the  same  direction  as  it  had  when  no  time  de- 
lay was  encountered  i.e.,  when  the  position  vector  was  0 P.  As  can  be  found 
this  occurs  when  q = ^ . Since  the  velocity  vector  in  the  systems  under  con- 

sideration  is  appreciably  greater  than  the  position  vector,  and  for  cases 

I 

where  the  trigger  point  is  only  slightly  shifted. 

As  the  analysis  proves  and  as  can  be  seen  by  comparing  the  three  plots 


WADCTR  53-202 


28 


presenting  the  velocity  curve  in  its  relative  position  to  the  trigger  point 
within  a whole  cycle  (Figures  4.4.1*  4*4*2,  and  4.4.3)  the  following  fact 
becomes  apparent: 

The  smaller  the  phase  angle  (in  radians)  between  the  zero  of  the 
velocity  component  and  the  trigger  point,  the  higher  the  frequency  of  the 
missile  oscillation.  Let  S be  this  phase  angle  (in  radians).  Then,  the 
relation  exists  ^ C'f'tij  C is  the  time  taken  for  the  surface  to  travel 

from  one  stop  to  the  other  stop. 

If  a noticeable  integrator  time  lag  is  introduced  into  the  vector 
plot  (position  vector  being  OP2),  the  time  line,  t*:t  - t*,  representing  the 
moment  when  triggering  occurs  for  this  system,  advances  relative  to  the 
triggering  time  line  for  a system  with  no  time  delay  at  all.  In  other  words: 
the  angle  C f becomes  smaller.  From  the  equation  above  cf*  can  become  smaller 
only  if  f increases.  The  vector  diagram.  Figure  4.4.5  is  drawn  for  a given 
frequency  (tt  = 8).  At  higher  frequencies  the  velocity  vector  will  have  a 
greater  magnitude  relative  to  the  true  position  vector  (OP).  The  result  of 
this  change  would  be  a slight  additional  advance  of  the  triggering  time  line 
which  would  in  turn  increase  the  frequency.  Although  this  appears  to  have  a 
cascaded  effect,  in  the  practical  case  the  phase  shifting  influence  of  the 
rate  gyro  opposes  such  a development  by  increasing  the  phase  lag  between  true 
angular  velocity  and  the  trigger  point;  this  prevents  the  system  from  a 
further  advance  in  this  direction  and  keeps  the  shift  of  the  triggering  time 
line  relatively  small. 

Consider  the  case  in  which  the  Time  Constant  is  just  large  enough 
that  the  position  vector  in  question  (let  it  be  Pm)  together  with  the  velocity 
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vector  (V)  forms  a summing  vector  (Rm)  which  has  the  same  direction  as  the 
vector  R in  a system  with  no  time  lag.  When  the  system  has  no  time  lag 
the  projection  of  the  position  vector  P on  to  the  trigger  time  line  t-t  is 
smaller  than  its  actual  length.  Since  the  time  line  rotates  clockwise  it 
can  be  seen  that  triggering  occurs  after  the  time  line  passed  the  maximum 
of  the  position  component.  However,  in  the  case  under  consideration,  the 
trigger  time  line  and  the  vector  Pm  are  in  coincidence  thus  indicating  that 
triggering  occurs  exactly  at  the  maximum  amplitude  of  the  position  component. 

Any  larger  time  delay  in  the  integrating  motor  component  connected  with  a 
further  advance  of  the  trigger  time  line  causes  the  system  to  trigger  before 
the  maximum  of  the  position  component  appears. 

With  increasing  values  of  the  time  delay  the  trigger  point  advances 
on  the  position  curve  from  a point  behind  the  maximum  through  the  peak  of  the 
position  component  to  points  ahead  of  the  maximum.  As  far  as  the  velocity  curve 
is  concerned  the  trigger  point  travels  only  on  the  ascending  linear  portion  of 
the  curve  ahead  of  a maximum:  point. 

The  foregoing  explanations  weie  made  by  assuming  the  velocity  compo- 
nent to  be  a sinusoid  and  hence  are  approximative.  In  an  ideal  symmetrical 
case  the  component  sensing  the  actual  missile  roll  velocity  produces  an  output 
signal  which  represents  the  angular  velocity  without  delay;  its  graphical 
picture  in  the  time  plot  is  composed  of  straight  lines  and  parts  of  a parabola 
hence,  the  deviation  from  a sinusoid  is  of  minor  magnitude.  In  addition,  if 
a component  simulating  the  actual  rate  gyro  performance  is  included,  the  output 
signal  then  assumes  a shape  which,  due  to  the  filter  effect  of  the  second  order 
component,  differs  from  a sinusoid  still  less. 
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The  vector  plot,  furthermore,  reveals  the  fact  that  the  frequency 
decreases  if  the  gain  of  the  position  component  is  increased.  The  vector 
V in  such  a case  has  to  be  reduced  (since  the  position  vector  is  assumed  to 
be  of  unity  length).  This  rotates  the  resulting  vector  R and  consequently 

r 

the  trigger  time  line  t-t  clockwise  yielding  a larger  angle  <*  and  therefore 
a smaller  value  for  the  frequency. 

For  very  large  time  lags  where  the  position  vector  almost  disappears, 
the  response  should  be  similar  to  a system  with  only  a velocity  signal  resulting 
in  an  almost  triangular  surface  motion  (see  general  remarks  in  introduction). 

The  deadspace  in  the  relay  sensing  is  jointly  dependent  on  the  fre- 
quency of  oscillation  and  the  magnitude  of  the  signal  (J,  “P  J this  is 

immediately  apparent  from  the  graphic  representation  of  this  deadspace  in 
Figure  4.4*6.  It  can  be  seen  from  careful  examination  of  the  mixer  signal 
in  the  vicinity  of  the  relay  deadspace  in  this  figure  that  the  delay  due  to 
relay  deadspace  is  a function  of  the  mixer  signal  (j,  V+ft  P)  slope  in  the 
deadspace.  This  is  in  turn  a function  of  the  frequency  and  amplitude  of  the 
mixer  signal  and  can  be  calculated  for  individual  cases. 

The  delay  in  relay  signal  due  to  relay  time  lag  is  not  frequency 
dependent.  For  the  relay  used  in  computer  analysis  it  was  approximately  0.02 
second. 

The  time  lag  in  the  actuator  is  also  not  rrequency  dependent  but  is 
characteristic  of  the  actuator  itself.  results  physically  in  a delay  in  the 
time  taken  by  the  surfaces  in  comcing  up  to  speed.  Once  a constant  speed  is 
reached  the  surfaces  continue  until  encountering  physically  stops  or  until  the 
relay  reverses. 


WADCTR  53-202 


31 


Auxiliary  lines  which  contribute  to  the  ease  of  application  have 
been  added  to  Figure  4.4.4.  In  the  velocity  curves,  lines  shewing  the  phase 
lag  caused  by  a rate  gyro,  ^ = 5 c.p.s.,  with  four  values  for  the  damping 
ratio,  are  plotted  against  the  various  frequencies.  It  is  characteristic  of 
the  plot  that  these  lines  are  nearly  straight  and  parallel  to  the  trigger  line. 
These  phase  lags  apply  to  related  position  curves  in  the  same  manner. 

Another  set  of  lines  has  been  added  in  the  position  half  of  the  plot 
to  indicate  phase  shifts  of  10°,  20°,  30°,  and  so  on  from  the  trigger  line 
for  the  various  frequencies.  These  will  be  used  with  Figure  4. 4.4*1  to  consider 
the  effects  of  time  lags  in  the  integrating  motor. 

The  utilization  of  Figure  4.4.4  can  be  best  demonstrated  by  con- 
sidering some  related  problems. 

Example  No.  1.  Consider  a system,  having  ideal  sensing,  in  which  because 
of  delay  in  the  discriminator  relay  0.02  sec  time  delay  occurs  in  the  triggering 
mechanism.  It  is  desired  to  maintain  the  system  at  1 c.p.s.  What  gain  ratio 
is  necessary? 

(1)  The  delay  results  in  a relative  shift  of  velocity  and  position 
curves  with  respect  to  the  acceleration  curve  and  hence  the 
triggering  line.  This  may  be  represented  in  Figure  4.4.4  as  a 
shift  of  the  trigger  line  equal  to  0.02  seconds  to  the  left. 

(2)  Then  reading  the  ratio  of  velocity  to  position  signals  on  the 
1 c.p.s.  curve  yields  velocity/position  equals  1.3/0.49  or 
1:0.377. 

(3)  The  necessary  gain  ratio  of  velocity/position,  to  be  adjusted 
in  a physical  setup,  therefore,  equals  1:2.65. 
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Example  No.  2.  Consider  a system  with  the  same  time  delay  in  the  relay 
and  a rate  gyro  with  fn  ■ 5 c.p.s.  and  ^ = 0.3;  which  is  to  oscillate  at 
0.7  c.p.s. 

(1)  Phase  lag  by  the  rate  gyro  at  0.7  c.p.s.  is  seen  to  be  repre- 
sented by  the  line  labeled^  = 0.3. 

(2)  Time  delay  treated  as  in  Example  No.  1,  except  that  it  is  now 
added  to  the  delay  contributed  by  the  gyro. 

(3)  Attenuation  or  amplification  of  the  true  velocity  signal  by 
the  gyro  is  unimportant  in  that  it  results  in  proportional 
changes  for  both  position  and  rate  signal. 

(4)  Referring  to  the  0.7  c.p.s.  curves  and  reading  the  ratio  of 
velocity/position  under  the  shifted  trigger  line  reveals  a ratio 
of  approximately  3:1. 

(5)  Hence  the  gain  ratio  of  velocity/position  necessary  in  the  physical 
set-up  to  maintain  0.7  c.p.s.  under  these  conditions  is  1:3. 

Example  Ho.  3.  Consider  a system  with  .02  second  delay  in  the  discriminator 
as  before  but  this  system  is  preset  at  a gain  ratio  of  1:2.  At  what  frequency 
will  it  oscillate  if  the  rate  gyro  ^ - 0.3  and  « 5 c.p.s.  is  used? 

(1)  Treat  the  time  delay  again  as  a shift  of  the  trigger  line  to  the 

left,  and  as  the  time  shift  by  the  rate  gyro  is  approximately 
the  same  for  all  frequencies  treated  here,  add  the  shift  indi- 
cated for  a rate  gyro  having  ^ - 0.3. 

(2)  Looking  under  the  new  trigger  line  find  the  frequency  which  has 
the  desired  ratio  at  triggering  in  this  case  approximately  0.95 
c + p*  s • 
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Example  No.  4.  Analogue  simulator  results  reveal  that  increasing  the 
time  delay  in  the  integrating  motor  does  not  have  a detrimental  effect  on 
the  steady  state.  Explanation  for  this  property  was  given  earlier  in  this 
section.  The  conclusion  drawn  can  be  substantiated  by  using  the  VPR  chart 
in  Figure  4.4.4.  Consider  a system  operating  with  the  following  characteris- 
tics: 

(1)  Angular  roll  rate  is  sensed  by  a rate  gyro  ^ = 5 c.p.s.  and 

j = 0.3. 

(2)  The  gain  ratio  of  rate  to  displacement  is  1:1. 

(3)  The  time  taken  for  the  Discriminator  Relay  signal  to  change 
sign  from  an  extreme  position  is  0.02  sec.  (Time  delay  in  the 
relay) . 

(4)  The  time  delay  in  the  integrating  motor  is  0.1  sec. 

Steps  in  determining  the  system  characteristic  oscillation  are  as  follows: 

(1)  The  effect  of  the  rate  gyro  is  to  move  the  ideal  tripping  line 
into  the  position  of  the  line  ^ m 0.3. 

(2)  A 0.02  sec.  delay  in  the  relay  results  in  a further  trigger  line 
shift  of  0.02  sec.  (see  time  scale)  to  the  left. 

(3)  By  comparing  the  frequency-related  pairs  of  lines  at  this  new 
trigger  line  it  is  found  that,  with  no  time  delay  in  the  inte- 
grator and  for  a gain  ratio  1:1,  the  system  would  oscillate  at 
1.1  c.p.s. 

(4)  Reference  to  the  associate  time  delay  chart  (Figure  4. 4. 4.1) 
reveals  that  a time  delay  of  0.1  sec  at  1.1  c.p.s.  results  in  a 
phase  shift  of  35  degrees  and  an  attenuation  of  0.83  for  the 
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position  signal 


(5)  The  new  tripping  line  maintains  its  position  in  the  velocity- 
plane  but,  since  the  time  delay  is  in  the  integrator,  moves 

an  additional  35  degrees  to  the  left  in  the  displacement  plane. 
Now  comparing  reading  at  the  new  tripping  point  yields  an 
attenuated  displacement  of  0.32°  and  a value  of  0.42°/sed, 
for  the  velocity  component  the  gain  ratio  thus  being  1:1.3 
at  this  point. 

(6)  However,  if  the  ratio  at  tripping  is  to  remain  at  1:1,  a rule 

of  thumb  established  for  this  system  can  be  applied:  the  lower 

the  ratio  the  higher  the  frequency.  Hence,  the  ratio  being 
lower  than  1:1.3,  a frequency  of  more  than  1.1  has  to  be  chosen. 

(7)  Resorting  now  to  the  method  of  iteration  yields  for  each  new 
frequency  considered  a new  phase  and  attenuation  due  to  the 
integrator  time  lag.  The  rule  of  thumb  in  paragraph  (6)  above 
eliminates  most  of  the  guess  work  from  this  process. 

(8)  Each  new  frequency  also  means  a new  phase  lag  due  to  the  rate 

gyro  effects.  However,  since  the  lines  referring  to  the  gyro 

effect  on  the  VPR  chart  are  represented  by  straight  line,  nearly 
parallel  to  the  trigger  line,  a small  change  in  frequency  re- 
sults in  only  a negligible  shift  of  the  trigger  line. 

(9)  At  an  estimated  freqxiency  of  1.15  c.p.s.  the  phase  lag  is  36 
degrees  and  the  attentuation  0.81.  Apparent  tripping  displace- 
ment at  this  point  is  (0.34  degrees  X 0.81)  0,28  degrees  as 
compared  to  an  apparent  tripping  velocity  of  approximately  0.25 
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degrees  yielding  an  apparent  gain  ratio  of  1:0.9,  this  being 
now  lower  than  the  previously  found  ratio  of  1:1,  indicating 
that  the  correct  frequency  will  be  slightly  below  1.15  c.p.s. 

An  estimate  of  1.14  c.p.s.  would  be  very  close;  increased 
accuracy  could  be  obtained  by  extending  the  scales  of  the  VPR 
chart.  Thus  the  introduction  of  this  integrator  time  lag 
resulted  in  an  increase  in  frequency  from  1.1  to  1.14  c.p.s. 

(10)  The  analogue  simulator  results  for  this  situation  yield  a fre- 
quency of  approximately  1.13  c.p.s.  Allowance  which  must  be 
made  for  VPR  chart  accuracy  limitations  are  such  that  this  agree- 
ment is  as  good  as  it  might  be  anticipated.  For  most  practical 
purposes,  this  degpee  of  accuracy  is  sufficient. 

The  VPR  chart  is  used  here  as  a graph  with  a supporting  graph  to  re- 
present attenuation  and  phase  lag  by  the  integrating  motor.  However,  this  chart 
would  be  much  more  convenient  for  practical  use  if  it  were  desired  as  a com- 
puter similar  in  appearance  to  the  E-4  navigation  computer  or  to  a slide  rule. 
Such  an  instrument,  applicable  to  a specific  system  could  be  designed  quite 
easily  and,  while  generalization  would  present  some  difficulty,  it  could  be 
accomplished  at  a comparatively  small  cost. 

4.5  Up  to  this  point  only  the  steady  state  operation  of  the  system  has 
been  considered.  If  the  system  is  workable  it  must  also  maintain  a stable  con- 
dition when  subjected  to  outside  disturbances.  It  would  appear  that  the  system 
could  react  to  a disturbance  by  varying  the  relative  time  of  dwell  at  the 
surface  stops.  However,  to  determine  the  details  of  system  performance  in  the 
disturbed  condition  by  theoretical  means,  it  would  be  necessary  to  conduct 
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an  additional  long-time  investigation  on  the  subject.  This  has  been  approached 
only  superficially  in  order  to  have  a basis  to  obtain  a comparison  for  analogue 
results.  A good  similarity  existed  and  it  was  decided  that  this  study  would  be 
completed  by  the  faster  analogue  methods. 

Consider  first  a well  mannered  system  of  this  type  in  which  the  center 
line  about  which  it  oscillates,  is  being  displaced  by  some  undetermined  force. 
Assume  only  that  when  displaced  in  such  a manner  the  system  maintains  its  fre- 
quency and  is  returned  to  its  neutral  position.  Figure  4.5*1  shows  the  missile 
position  in  thelypothetical  situation.  On  the  same  figure  velocity  and  accele- 
ration curves  were  built  up  in  a stepwise  manner  from  the  fact  that  V E 
Jk  dt  + c and  that  at  the  tripping  point  velocity  is  equal  to  displacement. 
Examination  of  the  figure  shows  that  this  disturbance  would  then  be  counter- 
acted by  varying  the  relative  dwell  time. 

Disturbances  are  referred  to  as  percentages  of  full  surface  deflection. 
An  external  disturbance  is  called  100#  if  a 100#  steady  state  surface  deflection 
will  exactly  counterbalance  it.  It  can  be  seen  that  if  the  system  is  to  be 
stable  under  the  influence  of  a disturbance  applied  in  a stepwise  manner  the 
resultant  acceleration  (disturbance  plus  surface  deflection  expressed  as  accele- 
ration) must  equal  zero.  Hence  the  average  acceleration  delivered  by  the 
control  surface  must  be  equal  in  magnitude  and  in  the  opposite  direction  from 
that  due  to  the  disturbance.  Then  the  line  of  average  acceleration  caused  by 
surface  deflection  is  merely  biased  from  the  original  zero  by  a percentage  of 
maximum  surface  displacement  equal  to  the  percentage  disturbance  considered. 
Figure  4.5.2  assiur.es  V- jAclt  + C V~fM  + and  shows  the  rela- 

tion which  must  exist  in  a stable  system  subjected  to  a light  disturbance  of 
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constant  magnitude.  It  is  of  particular  importance  to  note  that  the  missile 
oscillates  past  its  old  zero  position  but  the  average  position  is  now  somewhat 
changed. 

Figure  4.5.3  shows  the  result  of  applying  any  disturbance  strong 
enough  so  that  the  surface  dwells  on  only  one  side.  It  is  of  interest  to 
note  that  in  this  condition  the  theoretical  velocity  signal  is  zero  at  the 
second  tripping  point  T2*  Because  of  the  triggering  condition  and  the  fact  that 
position  maximum  occurs  at  m = o it  must  be  concluded  that,  in  the  theore- 
tical case  the  missile  will  always  return  to  its  original  zero.  Practically 
speaking  the  velocity  signal  would  indlude  a phase  lag  and  its  value  would 
differ  somewhat  from  zero  at  the  tripping  point.  Since  this  point  is  not  far 
from  the  peak  of  the  position  curve,  it  must  be  surmised  that  the  bomb  will  no 
longer  return  to  its  original  zero  position.  The  external  disturbance  condition 
which  will  cause  the  surface  to  dwell  at  only  one  of  the  stops  is  determined 
by  the  following  parameters:  frequency  by  which  the  missile  oscillates,  surface 

stops,  surface  speed,  and  surface  effectiveness.  These  are  related  by  simple 
geometric  considerations . 


SECTION  IV 
SIMULATION  RESULTS 

5.1  The  system  first  simulated  was  one  in  which  the  actuator  output 
was  of  an  isosceles,  trapezoidal  wave  form.  It  is  assumed  that  the  actuator 
output  is  independent  of  any  load  placed  upon  it  in  the  form  of  an  external 
hinge  moment  and  that  the  surfaces  will  be  limited  to  i 20  degrees  deflection. 
Due  to  the  method  of  simulation  used  for  the  actuator  the  surface  dwelling  at 
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the  stops  is  represented  as  a high  frequency,  low  amplitude  vibration  whose 
mean  value  is  the  constant  output  representing  true  surface  dwelling.  A 
comparison  of  the  desired  actuator  travel  to  the  actuator  travel  taken  from 
an  oscillogram  recording  of  the  simulated  actuator  appears  in  Figure  5.1.1. 

As  has  been  mentioned  earlier  in  this  report,  acceleration  is  directly 
proportional  to  aileron  deflection.  This  being  the  case  the  angular  velocity 
is  obtained  by  integrating  the  acceleration  curve.  A comparison  between  theo- 
retical angular  velocity  and  angular  velocity  taken  from  oscillogram  recording 
of  the  simulation  is  shown  in  Figure  5.1.2.  Similarly  the  angular  displacement, 
obtained  by  integrating  angular  velocity  appears  as  Figure  5.1.3,  with  that 
taken  from  an  oscillogram  recording  of  simulation. 

5.2  The  setup  shown  in  Figure  5.2.1  was  used  to  evaluate  the  effects  of 
varying  the  surface  angular  velocity  on  the  frequency  and  amplitude  of  missile 
oscillation.  Since  this  case  was  investigated  under  the  assumption  of  an  ideal 
sensing  device  the  excellent  agreement  between  the  results  of  the  simulation  and 
the  theoretical  evaluation  was  anticipated.  The  results  are  presented  graphically 
in  Figure  5.2.2  showing  frequency  and  amplitude  of  oscillation  vs  aileron  speed. 
Change  in  aileron  effectiveness  (missile  roll  acceleration  (deg/sec-^)  per 
degree  aileron  deflection)  is  seen  in  the  theoretical  analysis  to  have  no  effect 
on  frequency  but  it  effects  amplitude  of  oscillation  directly. 

A rate  gyro  y#  - 5 c.p.s.  and  J s 0,316  was  next  introduced  into  the 
system  as  the  means  of  obtaining  angular  roll  rate  of  the  missile.  The  setup 
used  for  the  following  investingation  is  shown  in  Figure  5.2.3.  The  effect  of 
varying  the  gain  ratio  under  this  condition  is  shown  in  Figure  5.2.4  which 
yields  frequency  and  amplitude  of  missile  oscillation  and  dwelling  time  of  the 
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ailerons  at  their  stops.  It  can  be  seen  that  increasing  the  gain  of  the 
positional  signal  results  in  a decrease  in  frequency  as  in  the  ideal  system 
considered  theoretically  in  paragraph  4*3j  the  mode  of  the  change  is,  however, 
somewhat  altered. 

This  setup  was  also  used  to  check  the  system  reaction  to  the  appli- 
cation of  an  external  disturbance.  Figure  5*2.5  reveals  the  upper  and  lower 
oscillation  peaks  of  the  system,  operating  at  3 different  gain  ratios,  to 
varying  amounts  of  external  disturbance  expressed  as  percentages.  Obviously 
disturbances  of  100$  will  result  in  complete  system  failure  and  hence,  no 
accelerations  in  excess  of  90$  were  considered.  The  frequency  of  oscillation 
as  a function  of  external  disturbance  is  shown  in  Figure  5*2.6  for  the  same 
three  different  gain  ratios.  Figure  5*2.6  reveals  that  the  frequency  reaches 
a partial  minimum,  increases  and  finally  drops  off  toward  zero  as  the  distur- 
bance is  increased.  Sven  in  these  extreme  conditions,  increasing  the  gain 
factor  of  the  positional  signal  decreases  frequency  (Figure  5*2.6). 

It  may  be  noted  that  the  displacement  curve  differs  more  and  more 
from  a sinusoid  as  the  size  of  the  external  disturbance  is  increased  as  shown 
in  the  recordings  in  Figure  5. 2. 6.1  and  5*2. 6. 2 for  63$  and  $0$  disturbance 
respectively.  This  is  to  be  expected  since  the  symmetry  of  the  acceleration 
is  progressively  more  distorted  as  the  disturbance  increases.  First  the  dwell 
time  at  the  stops  is  no  longer  equal,  then  the  surfaces  dwell  at  only  one 
stop  but  travel  past  the  neutral  position,  and  finally,  as  the  disturbance 
becomes  quite  large,  the  surface  no  longer  returns  to  the  opposite  side  cf 
neutral  but  completes  its  entire  oscillation  deflected  on  the  one  side  of 
neutral  which  is  in  opposition  to  the  disturbance.  It  is  possible  to  compute 
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the  mean  value  of  surface  deflection  for  each  test  with  various  parameters 
by  evaluating  the  records.  For  a disturbance  of  approximately  50%,  one 
would  expect  the  mean  value  of  surface  deflection  to  be  10°  (»  50%  of  total 
surface  deflection).  Figure  5.2.7  yields  the  curve  for  computed  mean  deflection 
compared  to  expected  mean  deflection.  However,  determination  of  the  mean 
deflection  from  the  computed  records  is  difficult  since  in  all  cases  the  width 
of  the  trace  itself  is  quite  significant  in  the  final  results.  Hence,  the 
deviation  between  expected  and  computed  mean  deflection  in  Figure  5.2.7  is, 
in  all  probability,  one  in  either  evaluation  or  the  recording  mechanism. 

This  viewpoint  is  substantiated  by  the  fact  that  there  is  no  apparent  shift 
in  the  angular  velocity  signal;  a drift  in  this  signal  is  inherently  related 
to  a realistic  deviation  of  the  type  shown.  Since  the  error  is  consistently 
in  one  direction  and  appears  only  when  the  surface  travel  becomes  asymmetrical 
to  the  point  that  the  surface  dwells  on  only  one  side  it  appears  to  be  an  error 
in  recording  rather  than  evaluation.  This  error,  probably  a time  lag  in  the 
pen  motor  of  the  Recorder,  apparently  cancels  out  in  more  symmetrical  record- 
ings but  in  asymmetrical  c ases  yields  a larger  apparent  dwell  time  than  the 
system  actually  has. 

The  effect  of  varying  the  damping  ratio  in  the  rate  gyro,  on  overall 
system  performance,  was  investigated  for  undisturbed  and  disturbed  flight 
conditions;  a 25%  acceleration  was  used  as  representative  of  a disturbance. 
Damping  ratios  were  varied  from  0.1  - 1.6;  the  results  are  plotted  in  Figure 
5.2.8  showing  amplitude  of  upper  and  lower  peaks  of  missile  oscillation  in 
terms  of  amplitude/aileron  effectiveness  and  also  showing  frequency  plotted 
against  damping  ratio  for  the  undisturbed  case.  Figure  5.2.9  presents  the  same 
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data  for  the  system  reacting  to  a disturbing  acceleration. 

5.3  It  has  been  revealed  by  previous  analysis  (4.4)  that  a time  delay- 
in  the  integrating  motor  should  have  a very  small  effect  on  the  overall  system 
performance.  This  conclusion  is  supported  by  the  simulation  represented  in 
Figure  5.3.1.  Results  of  the  simulation  appear  in  Figures  5.3.2  and  5.3.3. 

5.3.2  shows  maximum  and  minimum  peaks  plotted  against  percentage  disturbance 
for  q (q  = l/T)  of  1,  2,  5,  and  10  for  a gain  ratio  of  1:1.  Figure  5.3.3 
shows  the  same  plot  for  q ■ 1,5,  and  15  if  the  gain  ratio  is  1:0.5.  Figure 
5. 3. 2.1  and  5. 3. 2. 2 are  recordings  for  q g 1 and  q = 10  and  for  gain  ratio  1:1. 
Superimposed  on  the  trace  for  external  disturbance  is  a continuous  sequence  of 
time  marks  indicating  periods  of  1 sec  each.  These  records  reveal  that  no 
great  change  in  system  performance  results  from  the  introduction  of  even  re- 
latively large  time  delays  in  the  integrating  motor.  The  magnitude  of  change 
which  can  be  tolerated  is  dictated  by  the  desired  accuracy  and  hence  must  be 
determined  for  individual  missile  control  system  combinations. 

A time  delay  in  the  actuator  has  a greater  effect  on  system  perfor- 
mance. The  computer  setup  showing  how  this  delay  was  simulated  is  presented 
in  Figure  5.3.4.  Time  delays  of  0,1  second  and  0.05  second  were  introduced 
into  the  actuator  and  the  results  are  plotted  in  Figure  5*3*5  which  represents 
the  results  for  q r 10  (T  ■ 0.1)  and  q = 20  (T  ■ 0.05  sec)  and  shows  maximum 
and  minimum  peak  values,  and  frequency  vs  percentage  disturbance.  A comparison 
of  Figure  5.3.5  with  5.2.5  reveals  that  a time  delay  of  0.1  sec  introduced  into 
the  actuator  has  a very  noticeable  effect  on  system  performance.  For  a time 
delay  of  0.1  sec  and  a gain  ratio  of  1;05,  the  maximum  mean  deviation  from  the 
original  missile  neutral  position  is  10  times  the  aileron  effectiveness  factor 
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as  compared  with  2.4  times  the  aileron  effectiveness  factor  for  no  time 
delay.  The  peak  deviation  in  the  first  case  occurs  at  a 25-30%  disturbance 
and  in  the  latter  at  63%  disturbance.  Hence  if  a system  is  expected  to 
operate  through  the  entire  disturbance  range  such  a time  delay  is  very  de- 
trimental to  system  performance;  however,  for  limited  disturbance  conditions 
the  time  delay  has  much  less  influence. 

5.4  The  setup  in  Figure  5«4»1  was  used  to  determine  the  effects  of  a 
dead  zone  in  the  sensitivity  of  the  flip-flop  relay.  A dead  zone  in  a sensing 
device  is  taken  to  be  a region  of  insensitivity  within  its  range  of  operation. 
The  dead  zone  for  the  relay  is  electrically  considered  as  a region  symmetrical 
about  the  zero  point,  which  has  a width  of  ilOVin  the  simulator  setups;  this 
would  correspond  to  a signal  representing  10  degrees/sec  angular  velocity. 

More  specifically  this  means  that  at  any  time,  the  mixer  signal  is  within  this 
dead  zone,  the  relay  senses  zero;  when  the  signal  changes  sigi,  the  flip-flop 
action  will  be  delayed  according  to  the  width  of  the  dead  zone  and  the  rate  at 
which  the  signal  is  changing  as  it  passes  through  zero. 

Since  the  dead  zone  in  the  relay  results  in  phase-shifting  the  tripping 
point  an  amount  determined  by  the  slope  of  the  mixer  signal,  and,  hence,  by  its 
magnitude,  the  response  of  a particular  system  is  dependent  on  system  gain  from 
the  rate  tyro  to  the  actuator.  Therefore,  the  term  amplitude  per  effectiveness 
has  no  meaning  as  used  in  previous  graphs;  for  Figures  5.4.2  and  5.4.3  the 
actual  amplitude  of  displacement  and  frequency  of  oscillation  vs  percent  dis- 
turbance is  shown  for  an  effectiveness  of  8 degrees/sec2/ degree  surface  deflec- 
tion. These  tests  also  included  integrator  time  delays  of  .05  and  .20  sec 
respectively.  Since  a time  delay  in  the  integrator  has  been  shown  to  have  only 
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a small  effect  on  system  response,  the  similarity  between  Figures  5«4.2  - 

and  5.4.3  was  anticipated.  These  fibres  also  show  that,  although  amplitude 
of  oscillation  and  displacement  become  relatively  large,  the  system  is  still 
well  natured  for  disturbances  of  up  to  40  percent*  The  recording  shown  in 
Figure  5*4»2.1  presents  the  manner  in  which  the  dead  space  influences  the 
system  operation.  It  would  appear  that  the  rehy  dead  space  might  have  a 
considerably  greater. influence  upon  the  transient  reponse  of  the  missiLe. 

The  nature  of  the  transient  resulting  from  a time  gust  would  differ  consi- 
derably from  that  resulting  from  the  step  input  used  here.  Hence,  conclusions 
concerning  the  transient  would  be  subject  for  further  study.  It  is  of  interest 
to  note  that  the  nature  of  the  transient  for  "step  like"  disturbances  besides 
being  dependent  on  the  disturbance  will  also  be  dependent  upon  the  position 
of  the  control  surfaces  at  the  time  of  its  application.  Since  the  control  sur- 
faces are  in  constant  oscillation,  the  results  of  a study  concerning  the  tran- 
sient would  be  limited  to  placing  upper  and  lower  limits  upon  the  system  output 
during  the  transient.  Figure  5. 4. 2. 2 through  5. 4*2. 4 show  a ^umber  of  tran- 
sient responses  for  the  symmetrical  system  operating  under  various  conditions. 

5.5  The  influence  of  component  characteristics  which  usually  cause  a 
phase  lag  have  been  discussed.  It  is  also  likely  that  certain  components  may 
have  a slightly  asymretrical  response  to  a symmetrical  forcing  function.  These 
non-symmetrical  tendencies  can  be  present  in  the  sensing  components  (rate  gyro 
and  integrating  motor)  and  in  the  controlling  or  actuating  components  (dis- 
criminator, actuator,  control  surfaces,  and  the  missile  itself). 

Consider  the  case  where  there  is  an  asymmetry  in  the  rate  gyro.  If 
the  missile  were  in  a symmetrical  steady  state  oscillation,  with  no 
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aerodynamical  feedback  from  the  controlling  elements,  the  signal  of  the 
rate  gyro  would  be  an  asymmetrical  oscillation  with  respect  to  its  original 
null  position.  In  other  words,  the  mean  angular  velocity  yielded  by  the  gyro 
output  would  have  a steady  state  error  determined  by  the  degree  of  asymmetry 
involved.  However,  for  the  closed  loop  the  system  is  to  minimize  an  error  in- 
dicated by  the  rate  gyro,  so  that  an  indicated  mean  steady  state  angular 
velocity  of  zero  results.  Consequently,  the  missile  must  have  a fixed  rate 
of  rotation  superimposed  on  its  oscillatory  roll  motion.  This  fixed  rate 
exactly  balances  the  asymmetry  of  the  gyro.  This  results  in  the  indicated  and 
time  angular  velocity  and  position  relationship  shown  in  Figure  5*5.1.  Thus, 
the  amount  of  this  type  of  error  which  can  be  tolerated  in  any  given  situation, 
is  merely  such  that:  rate  gyro  error  x time  * allowable  displacement  error. 

This  situation  has  been  investigated  by  a REA.C  simulation  which  supported  both 
qualitatively  and  quantitatively  the  conclusions  just  drawn.  The  results  of 
this  simulation  are  recorded  in  Figure  5. 5. 1.1.  The  system  for  which  the  re- 
cording was  made  had  the  following  characteristics:  surface  speed  120°/sec, 

surface  stops  at  ± 20°,  surface  effectiveness  14°/ sec^  per  degree  aileron  de- 
flection, rate  gyro  = 5 c.p.s.,  ^ = 0.316.  The  asymmetrical  output  of  the 
gyro  is  assumed  to  be  such  that  the  signal  is  amplified  by  1.3  when  it  is  of 
one  polarity  and  attenuated  by  0.7  when  the  signal  is  of  the  opposite  polarity. 
To  calculate  the  effect  of  any  asymmetry  in  the  rate  gyro,  one  must  merely 
determine  the  true  angular  velocity  which  results  in  an  indicated  angular 
velocity  whose  integral  over  one  cycle  is  zero.  Since  they  do  not  contribute 
to  the  asymmetry  to  the  indicated  angular  velocity,  the  ordinary  phase  lag  and 
attenuation  of  the  rate  gyro  system  have  no  significance  in  these  calculations. 
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Consequently,  by  assuring  a missile  angular  velocity  of  the  form 
" /4  SOI  tut/  —J)  ",  the  following  condition  exists: 


f ^ 
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Where:  D is  the  average  missile  angular  velocity 

A is  the  amplitude  of  missile  angular  velocity  oscillation. 

The  amplification  of  rate  gyro  signals  of  opposite  polarity  has  the 
ratio  l:a. 

When  this  is  reduced  it  yields: 
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from  which  D may  be  calculated  for  a specific  case. 

For  the  case  represented  by  Figure  5.5 #H>  A was  35°  and  a = 1.85 
which  when  substituted  into  equation  (33)  yields  D~6. 8-  This  indicates  a 
rate  of  missile  drift  of  about  6.8°/sec  as  compared  to  a simulated  drift  of 
?.2°/sec.  Each  type  of  asymmetry  would  require  special  consideration  to 
establish  its  effect  on  system  performance.  Generally  speaking,  asymmetry 
in  a rate  sensing  element  will  result  in  error  in  missile  angular  velocity 
while  asymmetry  in  position  sensing  devices  will  result  in  errors  in  missile 
angular  position. 

If  controlling  components  such  as  actuators,  control  surfaces,  sur- 
face stops,  discriminators,  etc,  are  asymmetrical  in  their  performance,  the 
result  of  the  asymmetry  will  be  apparent  at  the  sensing  element  outputs 
causing  the  system  to  take  corrective  action  in  such  a manner  that  for  the 
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steady  state  oscillation  the  average  indicated  missile  angular  velocity 
is  zero  (if  the  rate  sensing  element  itself  does  not  have  an  asymmetrical 
output).  Thus  the  mean  value  of  the  true  missile  angular  velocity  will 
also  be  zero.  In  terms  of  missile  motion,  the  final  result  will  be  a 
steady  position  error.  In  any  of  these  cases  the  introduction  of  an  asym- 
metry would  be  likely  to  result  in  a change  in  frequency  as  well.  However, 
since  it  was  only  desired  to  determine  the  steady  state  positional  changes 
brought  about  by  the  introduction  of  asymmetries  into  the  system,  the  changes 
in  frequency  and  amplitude  of  oscillation  were  not  investigated. 

6.1  AS  has  been  previously  mentioned  (2.3)  it  is  possible  to  improve  the 
performance  of  the  flicker  control  system  to  a certain  extent  by  designing  the 
surface  actuating  component  so  that  a feathering  operation  results.  For  the 
feathering  type  the  outbound  surface  motion  occurs,  as  before,  with  constant 
angular  velocity.  The  inbound  motion,  however,  is  independent  of  the  systems 
performance  inasmuch  as  the  surface  is  released  from  the  driving  parts  of  the 
servo  actuator  as  soon  as  the  summing  command  changes  its  sigi.  At  this  moment 
the  surface  swings  back  into  the  neutral  (streamlined)  position  due  to  the  force 
of  the  airstream.  From  this  point  it  is  motor  driven  again  into  the  outbound 
direction  until  it  hits  the  stop  or  is  released  once  more.  Details  relating 
to  the  specific  kind  of  motion  and  the  manner  by  which  it  has  been  simulated 
on  the  REAC  computing  machine,  have  been  explored  in  3.2B  and  by  Figure  3*2.6. 

The  sensing  device  (rate  gyro)  used  in  the  Feathering  Case  is  the 
same  as  that  used  in  the  symmetrical  case.  However,  because  of  the  manner  in 
which  the  surface  is  returned  to  neutral  in  the  feathering  case  it  is  an  in- 
herently faster  system  than  the  other.  This  results  in  system  frequencies 
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nearer  the  natural  frequency  of  the  rate  gyro.  As  can  be  seen  from  a parti- 
cular  feathering  performance  graphically  represented  in  Figure  4.3«4>  the 
angular  velocity  curve  for  the  missile  differs  considerably  from  a sinusoid. 
However,  because  of  filtering  action  of  the  gyro,  the  indicated  angular  ve- 
locity signal  which  appears  at  the  mixer  and  the  integrator  is  more  nearly 
sinusoidal.  Since  the  mixer  signal  is  a major  component  in  the  determination 
of  syster  performance,  this  filtering  effect  becomes  very  significant  in  the 
overall  picture. 

It  is  necessary  to  determ  ine  the  manner  in  which  the  system  performance 
is  altered  by  the  introduction  of  a rate  gyro  as  a sensing  device;  particu- 
larly the  phase  lag  and  distortion  of  the  s.igials  entering  the  integrator  and 
mixer  will  be  significant  since  they  change  the  apparent  ratio  of  velocitv 
signal  to  position  signal  at  the  time  the  system  triggers.  A very  good 
approximation  of  missile  motion  and  gyro  response  to  such  motion  maty  be 
gotten  by  Fourier  methods.  Fourier  representation  of  the  acceleration  of  the 
symmetrical  type  system,  shown  in  Figure  4.4.1  and  discussed  by  paragraph  4.3 
by  equating  P - 1,  is  given  by  the  following  approximation: 

(34)  y (t/  * "jjr^  » oi  4A  Lit  + 4*4  3*  4*4  iut.  ■+  2s~  5V -4*4  ) 

Where  o(  r (ST and  A is  the  physical  limit  to  aileron  travel  in  degrees  if  aileron 
effectivity  is  l°/sec^/degree  aileron  travel. 

Fourier  representation  of  the  acceleration  of  the  Feathering  System., 
shown  in  Figure  4.3.4  and  discussed  in  paragraph  4.3  by  equating  M = 0,  is 
given  by  the  following  approximation: 


WADCTR  53-202 


LB 


(35) 


f 


w 


A U ^ ( 2/  /L  — A,tT J un  cj  ^ 


irlut~u 

•b  ~q  3 /ta  AaZ.  2>u)t  — J (f-2,l  — 0&->  3/ilT ^ On  Sco'tJ 


Where  " A-  " is  the  fraction  of  a half  cycle  during  which  the  surfaces  travel 

outbound,  i.e. , r,  "™  . 

fi,  + Kj 

In  order  to  obtain  a clearer  picture  of  the  gyro  effect  on  the  two 
systems,  a comparison  of  missile  motion  to  motion  indicated  by  the  sensing 
device  is  made  in  Figure  6.1.1,  6.1.2,  6.1.3,  and  6.1.4.  In  figure  6.1.1 
and  6.1.2  missile  angular  velocity  and  position  are  plotted  on  the  same 
graph  with  indicated  angular  velocity  and  position  for  the  symmetrical  system; 
at  CJ  — 3 and  a = 6.  In  Figure  6.1.3  and  6.1.4,  a similar  comparison  is  made 
for  the  feathering  system  at  Cj  = 6.28  and  oi  = 12.56. 

The  symmetrical  system  shown  in  6.1.1  and  6.1.2  was  based  upon  the 
following  consideration: 

1.  Physical  surface  stops  at  ^20° 

2.  Surface  velocity  between  stops  120°  per  sec 

3.  Surface  effectiveness  1°/ sec^/deg  surface  travel 

4.  The  rate  gyro  utilized  had  = 5 c.p.s.  and ^ = O.316. 

5.  The  integrating  motor  operated  ideally. 

Such  a system  has  a theoretical  maximum  frequency  of  1.5  c;p.s.  and  it  is 
slowed  down  by  the  phase  lag  in  t>  e rate  gyro.  True  missile  r otion  was  gotten 
from  the  formulae  in  4.3.  Indicated  missile  velocity  was  gotten  by  integrating 
termwise  the  Fourier  equation  (No.  34)  and  applying  the  dynamics  of  the  rate 
gyro  to  the  results  of  this  integration  by  means  of  frequency  response 
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charts.  Indicated  angular  position  is  then  obtained  by  termwise  integration 
of  indicated  angular  velocity. 

The  ratio  of  position  to  velocity  at  the  triggering  point,  taken 
from  6.1.1  and  6.1.2  for  W = 3*  and  w r 6,  is  1:5  and  1:3*9  respectively  for 
time  missile  motion  and  is  1:4.6  and  1:3.2  respectively  for  indicated  missile 
motion.  This  reflects  a decrease  in  the  magnitude  of  velocity  relative  to 
position  when  the  gyro  is  introduced;  however,  when  frequency  is  increased, 
the  ratio  of  position  to  velocity  at  the  tripping  point  decreases  for  both 
true  and  indicated  missile  motions.  Hence,  approximations  made  by  the  V.P.R. 
chart  in  Figure  4.4.4  yield  qualitatively  correct  results.  Further  examination 
shows  that  the  slopes  of  both  indicated  and  actual  missile  angular  velocity  in 
the  region  of  the  triggering  point  are  nearly  the  same.  Hence,  the  method  of 
delaying  the  actual  missile  motion,  relative  to  the  triggering  line,  to  obtain 
the  indicated  missile  motion  gives  reasonably  accurate  results. 

Similar  examination  was  made  of  the  feathering  system  in  Figure  6.1.3 
and  6.1.4  for  w = 6.28  and  M = 12.56  respectively.  Ratios  of  position  to 
velocity  were  taken  at  the  triggering  point  and  found  to  be  1:4.5  and  1:9 
respectively  for  actual  missile  motion  and  1:2.9  and  1:2.3  respectively  for 
indicated  missile  motion.  Again  the  introduction  of  the  rate  gyro  results  in 
a decrease  of  velocity  signal  relative  to  position  signal  at  the  triggering 
point.  However,  in  this  system,  as  the  frequency  increases  the  ratio  of  missile 
true  velocity  to  position  increases  at  the  triggering  point  while  d.he  ratio  of 
indicated  missile  velocity  to  indicated  position  decreases.  These  frequencies 
were  intentionally  held  below  those  at  which  the  feathering  system  tends  to 
operate  and  hence  present  the  best  gyro  representation  likely  for  this  case. 
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* Comparing  these  missile  angular  velocity  curves  to  indicated  angular  velocity 
in  the  triggering  region  reveals  that  no  amount  of  shifting  of  the  true  curve 
relative  to  the  trigger  line  could  render  it  a passable  approximation  for  the 
indicated  curve.  This  is  true,  basically,  for  two  reasons:  (l)  because  the 
feathering  system  operates  at  a higher  frequency  and  thus  nearer  the  natural 
frequency  of  the  rate  gyro  than  the  hon-feathering  and  (2)  because  of  lie  lack 
of  symmetry  in  missile  motions  for  this  system,  the  harmonics  play  a more 
significant  part  in  gyro  outputs.  Hence,  a chart  utilizing  theoretically  true 
missile  motion  as  is  done  by  the  VPR  chart,  is  not  feasible.  However,  a work- 
able chart  could  be  prepared  for  a given  gyro  by  using  Fourier  approximations 
as  they  have  been  used  here,  to  obtain  theoretical  indication  missile  motion. 
Preparation  of  such  a chart  would  obviously  have  to  be  made  in  conjunction  with 
a particular  system. 

6.2  Tests  have  been  conducted  simulating  a missile  equipped  with  a rate 
gyro  and  a surface  actuator  with  a dynamic  performance  as  described  above. 

The  setup  on  the  REAC  for  simulating  a missile  with  a feathering 
type  control  system  is  given  in  Figure  6.2.1.  As  can  be  seen,  the  general 
outline  (except  the  surface  motion  component)  is  the  same  as  was  used  in  the 
so  called  symmetrical  (non-feathering)  case.  This  setup  represents  the  simu- 
lation of  the  missile  performance  involving  no  components  which  cause  time 
delays  or  dead  zones  or  other  similar  influences.  The  setup  consists  of 
practically  the  same  circuitry  used  for  the  symmetrical  system,  with  the 
feathering  component,  described  by  Figure  3.2.6,  substituted  for  the  inbound 
surface  actuator  component. 

It  seemed  necessary  to  extend  the  investigation  in  such  a manner  that 
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the  influence  of  a time  delay  in  the  integrating  motor  or  in  the  actuator 
could  be  studied.  The  same  subcomponaits  (representing  the  simulation  of  such 
an  influence)  were  used  and  placed  within  the  general  circ’iit  diagram  as  has 
been  done  in  the  symmetrical  (non-feathering)  case.  Therefore,  no  specific 
circuit  diagrams  are  presented  as  their  arrangement  can  easily  be  composed 
from  the  diagrams  for  the  corresponding  cases  relating  to  the  non-feathering 
type. 

In  order  to  obtain  information  on  the  missile  response  and  the 
performance  of  certain  individual  components  in  close  correlation  with  the 
developmental  progress  of  the  missile,  all  tests  were  run  with  an  outbound 
surface  angular  velocity  of  120°/sec,  an  inbound  return  following  the  equation 
as  given  in  3.2.B,  and  a surface  effectiveness  of  15°/sec^  missile  roll  ac- 
celeration for  1°  surface  deflection. 

6.3  One  apparent  characteristic  of  the  missile  "feathering"  flicker- 
type-control  combination  is  the  rather  high  frequency  of  oscillation  and  con- 
sequently a smaller  amplitude  of  displacement. 

Preliminary  inve station  revealed  that  the  frequency  is  markedly 
influenced  by  the  phase  lag  existing  between  true  roll  angular  velocity  and 
indicated  angular  velocity  as  yielded  by  the  output  signal  of  the  rate  gyro. 
Since,  for  frequency  ranges  below  the  natural  frequency  of  the  rate  gyro, 
the  phase  angle  becomes  larger  for  higgler  ratios  of  the  rate  gyro  damping  and 
since,  for  a given  rate  sensitive  device,  the  damping  ratio  is  a part  of  the 
system  which  can  be  easily  adjusted  within  a certain  range,  it  was  of  interest 
to  investigate  the  influence  of  the  damping  ratio  on  missile  response.  The 
rate  gyro  investigated  had  a natural  frequency  of  5 c.p.s.  the  same  as  in 
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previous  investigations  concerning  the  non-feathering  case. 

Figure  6.3.1  shows  curves  representing  the  frequency  of  the  missile 
roll  motion  plotted  against  the  gain  ratio  (velocity  component:  position 
component)  as  it  has  been  done  for  the  non-feathering  type.  There  are  four 
curves,  three  of  which  represent  the  frequency  dependency  on  gain  ratio  for 
different  damping  ratios.  The  fourth  curve  pertains  to  a test  in  which  the 
gyro  with  a relatively  high  damping  ratio  (J>  - 0.316)  was  used  in  con, bination 
with  a lead  network.  Figure  6.3.2.  The  output  of  this  combination  was  then 
utilized  in  the  same  manner  as  the  gyro  output  in  the  systems  without  the  lead 
netxrork.  The  purpose  of  this  lead  network  is  to  compensate  for  the  phase  lag 
of  the  rate  gyro  output  and  thus  to  shift  the  missile  frequency  closer  to  the 
natural  frequency  of  the  rate  gyro.  The  transfer  function  of  the  network  has 
been  found  to  be: 


’Ltl 
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Insertion  of  the  lead  network  results  in  lessening  the  dependency  of  missile 
frequency  upon  the  gain  ratio.  The  missile  frequency  which  is  relatively 
high  varies  only  aanall  amount  for  large  changes  in  gain  ratio. 

Corresponding  to  the  frequency  carves  are  those  for  the  amplitude 
of  the  missile  roll  oscillation.  Utilization  of  a rate  gyro  with  a natural 
frequency  of  5 c.p.s.  and  a damping  ratio  of  ^ - 0.316  in  combination  with 
a proper  lead  network  results  in  a large  reduction  of  the  amplitude  for  the 
steady  state  oscillation  and  a lesser  degree  of  dependency  on  the  gain  ratio 
(gain  of  velocity/gain  of  position)  as  shown  in  Figure  6.3.3. 

In  respect  to  design  work  connected  with  the  rate  gyro  it  is  of 
importance  to  obtain  information  about  the  magnitude  of  the  angular  velocity 


WADCTR  53-202 


53 


encountered  during  the  missile  performance.  Under  the  assumption  that  the 
missile  oscillation  is  of  an  almost  sinusoidal  nature  the  response  charac- 
teristic of  the  rate  gyro  gimbal  can  be  computed  if  both  the  missile  roll 
oscillation  and  the  transfer  function  of  the  rate  gyro  are  known.  Such 
an  assumption,  however,  is  justified  only  as  an  approximation;  reference 
to  Figures  6.1.3  and  6.1.4  shows  that  less  accurate  approximations  may  be 
also  made  by  the  formuhe  in  4*3.  The  graphical  representation  of  the  half 
cycle  of  the  roll  acceleration  (originated  by  the  surface  motion)  is  a 
rather  asymmetrical  and  distrted  triangle,  and  the  integral  of  it  yielding 
the  missile  roll  velocity  also  differs  somewhat  from  a true  sinusoid. 

The  evaluation  of  tests  simulating  the  missile  behavior  for  dif- 
ferent parameters  of  the  damping  ratio  in  the  rate  gyro  component  is  repre- 
sented in  Figure  6.3.4;  the  amplitude  curves  for  true  angular  velocity  and 
angular  velocity  indicated  by  the  rate  gyro  are  plotted  against  the  gain 
ratio. 

It  may  also  become  necessary  to  determine  peak  values  of  angular 
velocity.  These  peaks  may  be  approximated  by  assuming  the  velocity  wave  form 
to  be  sinusoidal  and  basing  calculations  on  this  assumption.  This  can  be 
seen  by  evaluating  the  records  to  determine  velocity  peak  amplitude  as  function 
of  frequency.  Figure  6.3.5  shows  three  plots  of  the  amplitude  ratio  of  a rate 
gyro  (rate  gyro  output/rate  gyro  input  or  in  other  words:  indicated  angular 

velocity/true  angular  velocity)  versus  frequency  ratio.  Parameters  are  again 
three  different  values  for  the  damping  ratio  of  the  rate  gyro.  The  curves 
represent  the  analytical  computation  for  the  amplitude  ratio  versus  frequency 
ratio  of  a rate  gyro  (second  order  system)  excited  by  a sinusoidal  forcing 
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function.  The  corresponding  values  taken  from  the  records  are  marked  by 
crosses.  There  is  a relatively  good  agreement  between  the  computed  response 
for  sinusoidal  oscillation  and  the  actual  values  as  yielded  by  REAC  simulations. 

In  order  to  obtain  a complete  picture  of  how  the  amplitude  of  roll 
angle,  true  and  indicated  roll  angular  velocity,  and  the  frequency  of  the 
missile  response  change  if  the  damping  ratio  in  the  rate  gyro  is  varied  over 
a large  range,  the  results  of  corresponding  tests  are  compiled  in  Figure 
6. 3. 6.1,  6. 3. 6. 2,  and  6. 3. 6. 3 respectively.  A specific  adjustment  for  the  gain 
ratio  (velocity/position)  of  1:1  has  been  chosen  for  this  investigation  after 
it  became  apparent  that  under  certain  conditions  this  ratio  (or  a similar  one 
in  this  neighborhood)  would  likely  produce  the  best  results  with  regard  to  the 
stabilization  problems  involved.  However,  the  qualitative  overall  picture  as 
revealed  by  these  figures  will  remain  basically  the  same  for  other  gain  ratios. 

Assuming  the  missile  oscillation  in  undisturbed  flight  may  be  approxi- 
mated by  a sinusoid,  it  could  be  expected  that  values  for  the  damping  ratio 
above  0.7  would  produce  a curve  for  indicated  roll  velocity  that  lies  below  the 
curve  for  the  true  angular  velocity  (Figure  6. 3. 6. 2);  such  a rate  gyro  system 
is  in  or  beyond  the  critically  damped  condition  and  no  amplitude  amplification 
takes  place.  On  the  other  hand,  the  phase  lag  (produced  by  the  second  order 
system  characteristic  of  the  rate  gyro)  increases  relatively  fast  with  fre- 
quency, thus  forcing  the  missile  control  system  to  oscillate  with  a lower 
frequency.  The  actual  attentuation  of  the  rate  gyro  output  is  insignificant 
since  the  system  is  sensitive  only  to  the  zeros  of  the  summing  command  which 
is  composed  of  both  the  indicated  angilar  velocity  signal  and  its  integral  thus 
eliminating  the  influence  of  absolute  magnitude  of  the  summing  command. 
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For  rather  small  values  of  the  damping  ratio  in  the  rate  gyro  a 
higher  frequency  can  be  anticipated  since  in  this  case  the  phase  lag  in- 
creases slowly  with  frequency  until  the  natural  frequency  of  the  rate  gyro 
is  approached.  However,  this  being  the  case,  the  amplification  of  the  out- 
put signal  is  comparatively  large. 

With  the  damping  ratio  continuously  increasing,  the  mutual  in- 
fluence between  phase  lag  increase  and  change  in  amplification  or  attenti  a- 
tion  with  decreasing  frequency  causes  the  curve  of  the  indicated  angular 
velocity  to  first  decrease,  pass  a minimum,  and  then  climb  up  and  cross  the  curve 
of  true  angular  velocity;  the  curve  for  the  frequency  steadily  descends  and  those 
for  true  angular  velocity  and  the  roll  position  ascend.  These  curves  could  be 
of  considerable  significance  in  design  work  with  this  type  of  system. 

To  study  the  missile  response  to  external  disturbances,  a step  roll 
acceleration  of  different  magnitudes  has  been  applied  to  the  system.  Again, 
the  amount  of  the  external  acceleration  is  given  in  percent  of  the  product  of 
surface  deflection  (in  degrees)  and  the  effectiveness  (degree  sec“^  angular 
acceleration  per  degree  surface  deflection).  The  results  as  obtained  by  a 
simulation  for  which  a rate  gyro  with  a natural  frequency  of  3.16  c.p.s. 
and  a damping  ratio  of  0.25  was  used,  are  presented  in  Figure  6. 3. 7.1,  6. 3. 7. 2, 
and  6. 3. 7. 3 for  gain  ratios  of  1:1,  1:2,  and  1:4  respectively.  No  time  lags 
or  dead  spaces  were  assumed  to  be  present. 

Markedly  better  results  were  obtained  by  utilizing  a rate  gyro  with 
a higher  natural  frequency  (5  c.p.s.)  in  which  the  damping  ratio  was  0.316. 

The  results  as  evaluated  from  simulation  tests  are  given  in  Figure  6. 3. 8.1, 

6. 3. 8. 2,  6. 3. 8. 3,  6. 3. 8. 4 and  6.3.8. 5 for  gain  ratios  1:05,  1:1,  1:2,  1:4,  and 
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1:8  respectively.  Recordings  for  ratios  1:1  and  1:4  are  shown  in  Figure 
6. 3.8. 2.1  and  6. 3. 8. 4.1  respectively.  Comparison  of  these  graphical  plots 
with  those  as  obtained  for  corresponding  conditions  in  the  non-feathering 
(symmetrical)  case  presented  in  Figure  5.2.5  exhibits  two  major  significant 
differences. 

1.  As  to  the  course  of  the  mean  value  of  displacement,  there  exists 
no  peak  in  the  curves  for  the  feathering  case.  Generally,  an  increase  of  ex- 
ternal acceleration  results  in  a larger  amount  of  missile  roll  deviation. 

Such  a peak,  however,  is  always  present  in  the  symmetrical  (non-feathering) 
case. 

2.  For  a given  parameter  (gain  ratio  of  velocity/position)  and  a 
value  of  external  acceleration  for  which  the  non-feathering  system  has  not 
yet  reached  the  peak  of  mean  displacement,  the  oscillation  amplitude  is  much 
smaller  in  the  feathering  case,  as  is  the  mean  deviation.  However,  since  the 
curves  for  the  non-feathering  system  pass  through  a maximum  whose  absolute 
value  and  location  relative  to  external  acceleration  vary  with  the  parameter, 
i.e.,  the  gain  ratio,  there  are  regions  of  large  external  disturbance  where 
the  non-feathering  system  yields  a smaller  mean  deviation  from  zero  position 
than  does  the  feathering  type. 

The  extent  to  which  the  missile  performance  is  ihfluenced  by  a time 
delay  introduced  in  the  integrating  motor  component,  has  been  explained  in 
paragraph  4.4.  There  it  became  evident  that  an  integrating  motor  with  a re- 
latively large  time  constant  has  no  detrimental  effect,  but  improves  perfor- 
mances inasmuch  as  the  mean  deviation  of  the  missile  in  response  to  external 
accelerations  becomes  smaller  than  in  the  case  of  an  ideal  integrating  motor. 
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The  simulation  dealing  with  a symmetrical  case  applying  different  values  for 
the  time  constant  yielded  results  which  were  in  agreement  with  considerations 
previously  made  on  a more  qualitative  basis.  This  can  be  seen  by  comparing 
the  Figure  5.2.5  with  5.3.2.  A similar  response  with  respect  to  binding  of 
the  missile  nearer  to  its  zero  position  was  observed  in  the  feathering  case. 
Results  obtained  from  simulation  tests  with  a time  delay  of  100  m sec  and 
for  different  gain  ratios  are  presented  in  Figure  6. 3. 9.1,  6. 3. 9. 2,  and  6. 3. 9*3 
for  gain  ratios  of  1:1,  1:2,  and  1:4  respectively  which  compare  favorably  with 
those  given  in  Figure  6. 3. 8. 2,  6. 3. 8. 3,  6. 3. 8. 4 respectively.  However,  it 
should  be  kept  in  mind  that  all  these  curves  represent  only  the  steady  state  of 
missile  oscillatory  deviation.  There  are  more  or  less  noticeable  differences 
in  both  the  frequency  by  which  the  missile  oscillates  and  its  transient  be- 
havior; a few  examples  are  demonstrated  in  Figure  6. 3. 9. 4,  6. 3.9.5  for  gain 
ratios  of  1:1  and  1:4  respectively. 

In  agreement  with  results  obtained  from  the  investigation  of  the  non- 
feathering type,  a time  delay  in  the  actuator  pushes  the  response  curve  of  the 
missile  into  the  region  of  larger  deviations;  however,  the  difference  in  system 
response  between  a time  delayed  actuator  and  an  ideal  one  is  not  as  big  as  it 
is  in  the  case  of  the  corresponding  investigation  on  the  non-feathering  type. 

An  example  of  this  shift  is  given  in  Figure  6.3.10.1,  6.3.10.2,  6.3.10.3  and 
6.3.10.4  for  the  time  delays  of  .01,  0.02,  0.05,  and  0.1;  a gain  ratio  of  1:4 
was  used  for  each  of  these  figures.  Recordings  representing  various  components 
of  a system  using  a time  delayed  actuator  are  shown  in  Figure  6.3.10.1.1  and 
6.3.10.2.1  for  gain  ratio  1:4  and  disturbances  of  40£  and  time  delays  of  0.1 
and  0.01  sec  respectively. 
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SECTION  V 
CONCLUSIONS 

The  Flicker  Type  Control  System  is  particularly  adapted  to  pilotless 
flight  vehicles.  It  is  quite  versatile  and  maintains  well  controlled  os- 
cillations with  components  upon  which  it  places  liberal  requirements.  Des- 
pite its  non-linear  features  it  resembles  linear  systems  in  many  respects. 

In  addition,  it  lends  itself  very  well  to  electronic  simulation.  On  a 
single  REAC  unit,  complete  simulation  (for  one  axis)  was  performed  with  a 
small  amount  of  supporting  equipment. 

The  faster  operating  feathering  method  of  returning  the  surface  to 
neutral  position  by  the  force  of  the  airstream  results  in  higher  system  fre- 
quencies and  consequently  smaller  amplitudes  than  the  symmetrical  system  which 
has  equal  velocities  for  surface  outbound  and  inbound  travel.  It  seems  rea- 
sonable that  this  feathering  system  would  yeld  better  overall  system  perfor- 
mance. However,  if  the  larger  amplitude  of  oscillation  appears  to  be  accep- 
table, the  simpler  symmetrical  system  will  actually  yield  better  results  in 
response  to  gust  conditions  which  are  large  with  respect  to  maximum  surface 
effectiveness.  (See  Figure  5.2.5,  6. 3. 8.1,  6.3. 8.2  and  6.3. 8.3) 

The  steady  state  operation  of  a flicker  type  system  may  be  approximated 
by  simplified  graphic  methods  which  permit  the  introduction  of  many  of  the 
practically  difficult  component  characteristrics.  (See  paragraph  4.4.  and 
Figure  4.4.4) 

The  frequency  of  oscillation  for  a system  in  which  the  surfaces  contact 
stops  depends  in  a linear  manner  upon  the  surface  angular  velocity  and  upon 
the  limitations  to  surface  travel.  The  amplitude  of  oscillation  for  this 
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system  depends  linearly  upon  surface  effectiveness.  (See  Figure  5.2.2) 

The  frequency  of  a system  whose  surfaces  do  not  contact  stops  is  not 
dependent  upon  surface  travel  velocity.  However,  the  amplitude  of  oscillation 
for  this  system  depends  in  a linear  manner  on  both  surface  speed  and  effec- 
tiveness. 

Controlling  the  ratio  of  position  sigial  gain  to  angular  velocity 
signal  gain  provides  an  effective  means  for  controlling  system  frequency  with- 
in its  limits  of  operation.  It  serves  as  a means  of  obtaining  desirable 
system  performance  in  response  to  external  disturbance.  (See  Figure  5.2./*.  and 

6.3.3) 

Although  the  system  performance  depends  on  rate  gyro  characteristics  in 
an  expectable  manner,  the  rate  gyro  performance  is  a point  of  particular  in- 
terest. There  exists  a range  of  rate  gyro  damping  ratios  for  each  missile  con- 
trol system  combination,  which  by  minimizing  gimbal  deflection  offers  a desira- 
ble gyro  response.  (See  Figure  6. 3. 6. 2) 

In  general,  time  delays  in  the  system  have  the  usual  deteriorating  effect. 
However,  a certain  range  of  time  delays  in  the  position  sensing  device  tends  to 
improve  system  performance.  (Compare  Figure  6. 3. 8. 3 with  6. 3. 9. 3,  and  Figure 
5.3.3  with  5.2.5) 

The  system  remains  well  mannered  in  spite  of  relatively  large  time  delays 
in  any  of  its  components.  However,  if  the  components  introduce  errors  of  an 
asymmetry,  the  result  is  a steady  state  error  in  either  angular  position  or 
velocity  depending  on  the  nature  of  the  asymmetry  and  the  component  in  which 
it  is  located.  (See  paragraph  5.5) 
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GENERALIZED  DIAGRAM  CF  A FLICKER  TYPE  SYSTEM 


FIGURE  2.2.1 


FLICKER-TYPE  SYSTEM  INVESTIGATED 


TIQURX  2.3.1 


DISTURBANCE 


TO  RECORDER 


DETAILED  BLOCK  DIAGRAM  Of  MISSILE 


81AC  DIAGRAM  Of  MISSILE  SDTOLAIIOI 


WADGTR  53-202 


fIGURE  3.2.1 

63 


To  Mo-rcn*. 


Cf>  Gtyi*x>  Output  if*  lWTM*J<rtm. 

['0%4U«^VS«£]  T)('pt+'&'p+\<*')  = -Gllf  Iwt><CKTtT>  T^ouC>^- 

J lANqiUUKRVt.UOCtTV^ 

[Un  rr&>UL^uss^Jk.fe%*^axr>3 


"D  - TUcrfc.  Gi>raLO  Qtwr&Ki.  'D^'pua^c.v.vavkit  * Outi*utS»c,mm.* 

^ ^ 5 ^ [s'fcc]  f5um -S^c.'5tu>vrf 

y - aO/<KAp»<v«%1R>crio  |.Q 

*-  - <<  fe] 

u>^-  N AcrottKa- Fo.wauojKicy  OpTwiu  VjKtx>K.vyn»v«>  SysTiwi 

P - 

<j,  * F^sjerrot*.  £Q 

if  « .|  Ct.)  FouaM^  Fov4c.-t»o*ni 

DWAIL1D  BLOCK  DIAGRAM  07  RAW  GIRO 


-a.*ry 


olcTO* 


aWD  OoT^vjfy 
>-f •*  / -, 

(j)  * 


o»VdcI=  K. 

o.c  *■& 


RXAC  DIAGRAM  07  RAW  GIRO  SIMULATION 


WADCm  53-202 


7IGURX  3.2.2 

64 


v I <4L 

ip  = ~p~'Lf  - lKi*t>\c.Krrtcte  <Rouu  K.*4<a*v»v 

| ^ u 

<~P  - SokAHAtV4C^  C©KA*V\*VKJT>  ■*- 


DETAILED  BLOCK  DIAGRAM  07  MOTOR  INTEGRATOR  AND  MIXER 


Y-Y*+g,JlfMt 


oais  HiUio  mocln  coMPomim  _ j_ 

POSITION  COMPONENT  <J  z 

REAC  DIAGRAM  0?  MOTOR  INTEGRATOR  AND  MIXER  SIMULATION 


HOURS  3.2.3 


WADCTR  53-202 


65 


Output  ^ 
To  M\^>S\Ub 


ACTUATOR  OUTPUT  FOR  FEATHERING  TYPE  SYSTEM 


3.2.7 


^-k 


T =T tMk  CeMVrK.M*r 


DIAGRAM  Of  A COMPOHEHT  PRODUCING  A TIMS  DELAY 


0 s A.mq,la  Op 
V "VSjerr fcor « okj  OfShkFT 


MOTOR  INTEGRATOR  WITH  ( ) AND  WITHOUT  ( ) TIKE  DELAY 


IIGUBE  3.2.8 


WAECTR  53-202 


WADGTR  53-202 


Z.-Zi-Zjn/L;  >Zj+ 

Z®  = <_7.a- 

BEAC  DIAGRAM  07  DEAD  ZONE  SIMULATIOH 
7IGUBX  3.2.9 

71 


AIRFRAME  CONTROL  BY  CONVENTIONAL  TYPE  AUTOPILOT 
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FIGURE  4.1.1 


RELATIONSHIP  OF  MISSILE  MOTION  COMPONENTS 
FOR  A SYSTEM  UTILIZING  A COMBINATION  OF  ANGULAR 
VELOCITY  AND  ANGULAR  POSITION 
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Figure  4.4*5 
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MISSILE  COMPONENTS  CF  A SYMMETRICAL  SYSTEM  OPPOSING  A 
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SIMULATOR  SET  REFRESH'] TI NG  SYMMETRICAL  SYSTEM  WITH  RATE  GYRO  SENSING 
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REAC  SIMULATION  01'  SYMMETRICAL 
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Percent  of  Pull  Surface  Deflection 
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Acceleration  in  percent  of  full  surface  deflection  times  effectiveness 


SYMMETRICAL  CASE 

Time  Lag  in  Integrating  Motor  (T  = 0.]  sec) 
Dead  Space  in  Sunning  Relay  ±10  Veits 
Ratio  1:0.5 

25%  External  Disturbance 
Surface  Effectiveness  8 degrees/sec2, 
per  degree  surface  deflection 


SYMMETRICAL  CASE 

Time  Lag  in  Integrating  Motor  (T  = 0.1  sec) 
Dead  Space  in  Summing  Relay  ± 10  Volts 
Ratio  1:0.5 

Surface  Effectiveness  8 degreas/sea^ 
per  degree  surface  deflection 
62.5  External  Disturbance 


Fig.  5. 4.2. 2 
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SYMMETRICAL  CASE 

Time  Lag  in  Actuator  (T  = 0.1  sec) 
Ratio  1:0.5 

25$  External  Disturbance 
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surface  effectiveness 
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SYMMETRICAL  CASE 
RATIO  1:0.5 

No  Time  Lag,  No  Dead  Zone 
25%  External 
Disturbance 
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FEATHERING  SYSTEM 
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RELATIVE  ± AMPLITUDE  OF  ROLL  OSCILLATION 
VERSUS  EXTERNAL  ACCELERATION  FOR  DIFFERENT 
COMMAND  RATIOS  (FEATHERING  SYSTEM) 


ACCELERATION  IN  PERCENT  OF  FULL  SURFACE  DEFLECTION 
( X • EFFECTIVENESS  = 100%) 


RELATIVE  ± AMPLITUDE  OP  ROLL  OSCILLATIOH 
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ACCELERATION  IN  PERCENT  OF  FULL  SURFACE  DEFLECTION  TIMES 


RELATIVE  l1  AMPLITUDE  OP  ROLL  OSCILLATION 
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ACCELERATION  IN  PERCENT  OP  PULL  SURFACE  DEFLECTION  TIMES  EFFECTIVENESS 


ACCELERATION  IN  PERCENT  07  PULL  SURFACE  DEFLECTION  TIMES  EFFECTIVENESS 


ACCELERATION  IN  PERCENT  OP  PULL  SURPACE  DEFLECTION  TIMES  EFFECTIVENESS 


RELATIVE  * AMPLITUDE  OF  ROLL  OSCILLATION 
VERSUS  EXTERNAL  DISTURBANCE 


WADCTR  53-202 


1U2 


ACCELERATION  IN  PERCENT  OF  FULL  SURFACE  DEFLECTION  TIMES  EFFECTIVENESS 
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ACCELERATION  IN  PERCENT  OP  PULL  SURFACE  DEFLECTION  TIMES  EFFECTIVENESS 
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acceleration  in  percent  of  full  surface  deflection  times  effectiveness 

FIG.  6.3.10.4 
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